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Abstract. We deal with quasistatic evolution problems in plasticity with softening, in 
the framework of small strain associative elastoplasticity. The presence of a nonconvex 
term due to the softening phenomenon requires a nontrivial extension of the variational 
framework for rate-independent problems to the case of a nonconvex energy functional. 
We argue that, in this case, the use of global minimizers in the corresponding incremental 
problems is not justified from the mechanical point of view. Thus, we analize a different 
selection criterion for the solutions of the quasistatic evolution problem, based on a vis- 
cous approximation. This leads to a generalized formulation in terms of Young measures, 
developed in the first part of the paper. In the second part we apply our approach to 
some concrete examples. 



Keywords: plasticity with softening, quasistatic evolution, rate-independent processes, shear 
bands, incremental problems, viscous approximation. Young measures. 

2000 Mathematics Subject Classification: 74C05 (74C10, 28A33, 74G65, 49J45, 35Q72) 



Contents 



1 . Introduction 

2. Notation and preliminary results 

2.1. Mathematical preliminaries 

2.2. Mechanical preliminaries 

3. The minimum problem 

4. Regularized evolution 

4.1. Definition and properties 

4.2. Proof of the existence theorem 

5. Approximable quasistatic evolution 

5.1. Definition and properties 

5.2. The existence result 

6. The case of simple shear 

7. A spatially homogeneous example 

7.1. The special form of the solution 

7.2. A second order equation for the regularized evolution 

7.3. Fast dynamics 

7.4. Complete description of the quasistatic evolution 

8. Examples with concentrations and oscillations 

8.1. Strain localization 

8.2. Oscillation of the internal variable 
References 



2 
5 
5 
8 
12 
14 
14 
17 
25 
25 
29 
31 
33 
33 
35 
38 
40 
49 
50 
53 
57 



Preprint SISSA 38/2006/M (June 2006) 

1 



2 



G. DAL MASO, A. DESIMONE, M.G. MORA, AND M. MORINI 



1. Introduction 

In plasticity theory the term softening refers to the reduction of the yield stress as plastic 
deformation proceeds. Classically this is described by a family of yield surfaces depending 
on a parameter C . The evolution laws are formulated in such a way that the yield surface 
shrinks when the time derivative of the plastic deformation is not zero (sec [14], [15], [17], 
and [19]). 

We deal with this problem in the quasistatic case, in the framework of small strain asso- 
ciative elastoplasticity in a bounded and Lipschitz domain il C K.'' , d > 2. For simplicity 
we consider only the case of no applied forces and of prescribed boundary displacements on a 
closed subset Fq of the boundary dfl with positive {d — 1) -dimensional measure. The lin- 
earized strain Eu , defined as the symmetric part of the spatial gradient of the displacement 
u, is decomposed as the sum Eu = e+p, whore e and p arc the elastic and plastic strains. 
The stress a is determined only by e , through the formula a — Ce , where C is the elasticity 
tensor. We assume that, for every value of the parameter the elastic domain - the set of 
admissible stresses enclosed by the yield surface - has the form {a G M^^^ : an € K{Q}, 
where Mf^^ is the space of symmetric dxd matrices, an denotes the dcviatoric part of cr, 
and if (C) is a subset of the subspace M^^'^ of trace- free symmetric matrices. To simplify 
the mathematics of the problem, we assume that the set 

K := {((7,0 e M^^'^xM : a G KiC)} 

is a compact convex neighbourhood of (0, 0) 

To express the evolution laws, it is convenient to introduce an internal variable z , related 
to ( by the equation ( = —V'{z), where R ^ R is a given function of class with 
bounded second derivatives, called the softening potential. 

The strong formulation of the quasistatic evolution problem consists in finding functions 
u{t,x), e{t,x), p{t,x), cr{t,x), z{t,x), and C{t,x) satisfying the following conditions for 
every t G [0, -|-oo) and every x G Cl: 

(sfl) additive decomposition: Eu{t, x) = e{t, x) + p{t, x) , 

(sf2) constitutive equations: (T{t, x) = Ce{t,x) and ^{t,x) = —V {z{t,x)) , 

(sf3) equilibrium: div<T(i, x) = 0, 

(sf4) stress constraint: cr{t,x)D G K{^{t,x)) , 

(sf5) associative flow rule: {p{t,x), z{t,x)) G NK{o'{t,x)D,C{t,x)) , 

where dots denote time derivatives and Nk{o',C) is the normal cone to K at (cr, C) in 
M^^'^xR. The evolution is driven by a prescribed time-dependent boundary condition 

u{t, x) = w{t, x) for every t G [0, +oo) and every a; G Fq . 

It is supplemented by initial conditions at t = and by the traction-free boundary condition 

cr(t, x)n{x) = for every t G [0, +oo) and every a; G Fi , 

where n{x) is the normal to dfl at x and Fi =: dfl \ Fq . 
Introducing the support function 

H{^,e):= sup i^-.a + eO, 

where the colon denotes the scalar product between matrices, the flow rule (sf5) can be 
written in the equivalent form 

(sf5') dissipation pseudo-potential formulation: {(T{t, x)d, C(*) x)) € dH{p{t, x), z{t, x)) , 

where dH{^, 6) denotes the subdiff'erential of H at 9) . 

If V is strictly convex, this model describes plasticity with hardening, where the yield 
surface expands when the time derivative of the plastic strain is not zero. In this case it 
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is possible to give a variational formulation of the problem, which does not assume the 
existence of time derivatives, and is based on the energies 

Q(e) := i / Ce{x):e{x)dx, n{p,z):^ / H{p{x), z{x)) dx , 
Jn Jn 

V{z) := / V{z{x))dx. 
Jn 

The term z) is used to introduce the notion of dissipation of a function t {p{t), z{t)) 
on an interval [a, b] C [0, +00) , defined by 

fc 

VH{p,z;a,b) := sup^n{p{tj) - p{t j-i), z{tj) - z{tj-i)) , 

where the supremum is taken over all finite sequences (tj) such that a = to < ti < ■ ■ ■ < 
tk-i <tk = b. 

According to the energetic approach to rate- independent processes developed in [24], 
[18], [23], the variational formulation of the quasistatic evolution problem consists in finding 
functions u{t,x), e{t,x), p(t,x), and z{t,x) satisfying the following conditions: 

(vfl) global stability: for every t G [0, +00) we have Eu{t) = e{t) + p{t) on Q, u{t) = 
w{t) on To , and 

Q(e(t)) + V(z(t)) < Q(e) + n{p - p{t), z - z{t)) + V{z) 

for every u, e, p, z such that Eii = e + p on O, u = w{t) on Fq; 
(vf2) energy inequality: for every T G [0, +00) we have 

Q{e{T))+VH{p,z;0,T) + V{z{T))<Q{e{0))+V{z{Oj)+ [ {a{t). Ew{t)) d.t . 

Jo 

In the convex case, thanks to the Euler conditions, (vfl) is equivalent to the following 

property: 

(vfl') stability: for every t G [0, +00) we have Eu(t) = e{t) +p{t) on fi, u{t) = w{t) on 
To, and 

diy<T{t,x) — for X e , cr(f, x)n(x) = for a; € Fi , 

{(T{t.x)D.C{t-x)) e K ioixen, 

where cr and ^ arc defined by (sf2). 

There is a vast literature on variational methods in the study of evolution problems in 
elasto-plasticity. Among the papers which are closest in spirit to our approach we quote 
[25], [26], [21], and [22]. 

In this paper we assume that V is concave, which reflects the fact that the yield surface 
shrinks as p{t, x) ^0. To simplify the mathematics of the problem, we also assume that the 
image of —V is contained in the interior of the projection of K onto the C-axis. By lack of 
convexity, condition (vfl') is no longer equivalent to (vfl), and we have only (vfl) =4> (vfl') . 
In this case the selection criterion provided by global minimality is not justified from the 
mechanical point of view. Indeed, as we shall show in [7], global minimality leads to missing 
the softening phenomenon altogether. 

We explore a different selection criterion, based on the approximation by solutions of some 
regularized evolution problems, depending on a small "viscosity" parameter e > 0. In its 
strong formulation this regularized problem consists in finding functions Ue{t,x), es{t,x), 
p^{t,x), (Te{t,x), Zs{t,x), and (^^{t,x) satisfying (sfl), (sf2), (sf3), and 

(sf4)£ regularized flow rule: {p^it, x), Ze{t, x)) = Nl^{a-e{t, x)d, Cdt, x)) , 
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where 




0) 



Pk being the projection onto K. li z and ^ were not present, this condition would 
coincide with the flow rule of Perzyna viscoplasticity. We observe that condition (sf4)£ is 
closely related to (sf5). It also acts as a penalization which leads to the stress constraint 
(sf4) as 0. 

The existence of a solution to the e -regularized evolution problem is proved by a vari- 
ational method based on time discretization and on the solution of suitable incremental 
minimum problems (sec Theorem 4.4). Some parts of the proof are inspired by [29]. The 
uniqueness is based on Gronwall's lemma. In Theorem 4.3 we also prove that this solution 
is characterized by the following conditions: 

(rel)^ equilibrium condition: for every t G [0, +00) we have Eus{t) = e^lt) + p^{t) in ft, 
Ue{t) = w{t) on Fq, and 

div (Tg{t,x) = ior X G , (Tg {t, x)n{x) = for x G Fi , 
(cTf (t. x)d - ep^{t, x), x) — ez^{t, x)) e K for x € fl , 

where (Te and Cs ^-re defined by ae{t,x) := Cee{t,x) and <^^{t,x) := —V'{ze{t,x)); 
(re2)e energy equality: for every T G [0, -|-oo) we have 



By accepting only those solutions of (vfl') and (vf2) which can be approximated by solu- 
tions of (sfl), (sf2), (sf3), {si4)^ (or, equivalently, by solutions of (rcl)^ and {rc2)^), we regard 
quasistatic evolution as the limiting case of a viscosity-driven dynamics (Definition 5.1). A 
similar approach in finite dimension was used in [9] . Other rate-independent problems with 
nonconvex energy have been studied in [3], [16], and [12]. The last paper considers a different 
regularizing term based on the space gradient of the internal variable. 

The main difficulty in our approach is due to the fact that, by the nonconvexity of the 
energy, the components and z^ of the solutions of the e -regularized problem may develop 
stronger and stronger space oscillations as £ ^ 0. As a consequence of this fact, their weak 
limits do not satisfy, in general, (vfl') and (vf2) (see Section 8). To overcome this difficulty, 
we propose a weaker formulation in terms of Young measures. 

Since the functional Ti and V have linear growth, the classical notion of Young measure 
is not enough. To take into account possible concentrations at infinity, we use the notion of 
generalized Young measure introduced in [8], [1], and [11], following the presentation of [6]. 

In addition, to write the Young measure version of (vf2) we need to introduce a notion of 
dissipation for a time-dependent family of generalized Young measures. A natural definition 
can be given by taking the limit of the dissipations of suitable time-dependent generating 
functions. Unfortunately, this limit does not depend only on the values of the generalized 
Young measures at each time, but it also involves the mutual correlations between oscillations 
at different times. We solve this problem by using the notion of system of generalized Young 
measures introduced in [6] . This allows us to write a Young measure formulation of problem 
(vfl'), (vf2) (Theorem 5.4) and to prove an existence result (Theorem 5.6). 

The second part of the paper is devoted to some examples. The first example, developed 
in Section 7, deals with a spatially homogeneous case, where the £-regularized evolution 
is described by a system of ordinary differential equations, and the study of the limit as 
£ ^ reduces to the analysis of a singular perturbation problem for ordinary differential 
equations. In this example the stress cr{t) converges, as t ^ cxd, to a constant matrix a^o , 
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which coincides with the yield stress that would be obtained in the perfectly plastic case 
with elastic domain 

In other words, as time tends to infinity, the material behaves in the weakest way permitted 
by its internal variable. Moreover, after a critical time, \cr{t)\ is decreasing with respect to 
t, reflecting the fact that the material softens as plastic deformation proceeds (see Figure 1). 
Finally, for some values of the parameters, this example exhibits a jump discontinuity at a 
certain time (see Figures 2 and 3). This leads to a strict inequality in the energy balance 
(vf2), which shows that an instantaneous dissipation occurs at the discontinuity time. 

A second group of examples, developed in Section 8, shows that strain localization, in 
the form of a shear band, may occur in this model even if the initial and boundary data 
are sufficiently regular. One of the examples exhibits also a strong oscillation of the internal 
variable localized near the shear band. As s ^ 0, this leads to a Young measure 
solution of the quasistatic evolution problem. In this example the usual weak* limit of Ze 
(always given by the barycentre of the Young measure solution) still satisfies the equilibrium 
condition (vfl'), but does not satisfy the energy inequality (vf2). This shows that, because 
of the lack of convexity, some important terms generated by the space oscillations of the 
approximate solutions can be captured only by the Young measure formulation. 

2. Notation and preliminary results 

2.1. Mathematical preliminaries. We begin with a quick presentation of the mathemat- 
ical tools used in the paper. 

Measures. The Lebesgue measure on M'', d > 1, is denoted by , and the A: -dimensional 
Hausdorff measure by H'^ . Given a Borel set S c M'' and a finite dimensional Hilbert 
space S, Mb(i3;S) denotes the space of bounded Borel measures on B with values in S, 
endowed with the norm := |/x|(iJ), where S Mb{B) := Mi,{B;R) is the variation 

of the measure /x. The space of nonnegative bounded Borel measures on B is denoted by 
M^{B). For every ji G Mb{B;E) we consider the Lebesgue decomposition /i = /i" + /i* , 
where /k" is absolutely continuous and /k* is singular with respect to Lebesgue measure . 

If /i* = 0, we always identify fi with its density with respect to Lebesgue measure 
In this way L^{B-,'E.) is regarded as a subspacc of Mi,(B:E), with the induced norm. In 
particular /j."- G {B; S) for every /j, G Mi,{B; S) . The norm, 1 < r < oo , is denoted by 
II • \\r- The brackets (•,•) denote the duality product between conjugate L'" spaces, as well 
as between other pairs of spaces, according to the context. The symbols V and A denote 
the maximum and minimum of two numbers or functions, while the symbol (•)"*" denotes 
the positive part. 

If B is locally compact (in the relative topology), by the Riesz representation theorem 
(see, e.g., [28, Theorem 6.19]) Mb{B;'E.) can be identified with the dual of Co(B;S), the 
space of continuous functions (p: B — > S such that {\ip\ > e} is compact for every £ > 0. 
The weak* topology of Mb{B; S) is defined using this duality. 

Matrices. The space of symm,etric dxd matrices is denoted by Mf^,J^ ; it is endowed 
with the euclidean scalar product C : C •= tr(^C) = ^ij^ijdj ^.nd with the corresponding 
euclidean norm j^j := {^-O^^^- The symmetrized tensor product a© 6 of two vectors 
a, b E is the symmetric matrix with entries (a,6j + ajbi)/2. It is easy to see that 
tr(a Qb) = a-b, the scalar product of a and b. 

When d > 2 we define M^'''' as the space of all matrices of Mf^^ with trace zero. It 
turns out that M^^** is the orthogonal complement of the subspace MJ spanned by the 
identity matrix I. For every ^ G Mf^^ the orthogonal projection of ^ on M.I is 2'tr{^)I , 
while the orthogonal projection on M^^'' is the deviator of ^, so that we have the 
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orthogonal decomposition 



The case d = 1 is special: we have M^^^ = K and we do not need any orthogonal decompo- 
sition of the space. For the purposes of this paper (sec Section 6) it is convenient to define 
M.]^^ := R and := ^ for every ^ G M, although this does not agree with the definition 
given for d>2. 

Functions with bounded deformation. Let U be an open set in R'' , d>l. For every 

u S L^{U; M.'^) let Eu be the Mf^,J^ -valued distribution on U , whose components are defined 
by EijU = {DjUi + DiUj)/2. The space BD{U) of functions with bounded deformation is 
the space of all u e V-{U\ W^) such that Eu e Mb{U; M^^^) . It is easy to see that BD{U) 
is a Banach space with the norm 



It is possible to prove that BD{U) is the dual of a normed space (see [20] and [31]). The 
weak* topology of BD{U) is defined using this duality. A sequence Uk converges to u 
weakly* in BD{U) if and only ii Uk ^ u weakly in L^{U;W^) and Euk Eu weakly* in 
M(,(?7; Mg^^) . Every bounded sequence in BD{U) has a weakly* convergent subsequence. 
Moreover, if U is bounded and has Lipschitz boundary, every bounded sequence in BD{U) 
has a subsequence which converges weakly in L'^/^'^-^\U;M.'^) and strongly in U{U\m.'^) for 
every r < d/ {d — 1) . For the general properties of BD{U) we refer to [30]. 

In our problem u G BD{U) represents the displacement of an elasto-plastic body and 
Eu is the corresponding linearized strain. 

Generalized Young measures. As mentioned in the introduction, the results of this paper 
are based on the notion of systems of generalized Young measures developed in [6] . For the 
reader's convenience we collect here the main definitions and our notational conventions, 
while we refer to [6] for the motivations behind them and for the main properties. 

Let [/ be a bounded open set in K'' , d>l. We shall apply the notions introduced in [6] 
with X := U and \ = . To define the space of generalized Young measures we introduce 
the space C''°'"(?7xS) of all / e C(f/xS) such that /(.x, •) is positively homogeneous of 
degree one on S for every x . This space is endowed with the norm 



where Eh := G S : |^| = 1}. The dual of the Banach space C'*°™(t/x5) is denoted by 
M*({7xS), and the corresponding dual norm by || • ||*; the weak* topology of M*({7xS) 
is defined by using this duality. As in the case of Af5(J7xS), it is sometimes convenient to 
write the dummy variables explicitly also for this duality product and to use the notation 
instead of 

The space SxR is endowed with the product Hilbert structure. The corresponding spaces 
C''°"([7x(SxIR)) and M,([7x(SxR)) are denoted by C'^'"(c7x5xR) and M,(C7xSxM). 
The space GY{U ; S) of generalized Young measures on U with values in S is defined as 
the set of all ji e M*(C/xSxlR) satisfying the following properties: 

(a) positivity property: 



||?i||i + ||£;w||i. 



WfWhom :=max{|/(a;,OI ■ x U , C e Sh} , 



(/, /x) > for every / e C''°'"(i7xSxM) with / > ; 



(2.1) 



(b) 



support property: 



(c) 



(/, At) = for every / e C^""" {U x'E.xB) vanishing on J7xSx [0, +oo) ; 
projection property: 




for every G C{U) . 



(2.2) 
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It is easy to see that for every /i e GY{U ; ^) we have 

M*^{VWTW,K^,tv))- (2.3) 

We recall that for any ^ G GY{U;E) the barycentre of fj., denoted by bar(/x), is defined 
as the element of Mb{U; S) such that 

bar(^)) = (^(x) . Kx, V)) (2.4) 

for every (f G C(L'' ; S) . \i ^ 1^ weakly* in GY{U; S) , then 

bar(/ifc) ^ bar(/i) weakly* in Mi,{U; S) . (2.5) 

Given a measure p G Mf,([/;S), the generalized Young measure associated with p is the 
element Sp of GF(C7;S) defined for every / e C''°'"(C7xSxM) by 

{f,Sp} := l_f{x,%{x),^{x))d\{x), (2.6) 

where A G Mj^{U) is an arbitrary measure such that « X and p << A (the homogene- 
ity of / implies that the integral does not depend on A). 

Given another finite dimensional Hilbert space S', let V' Ux'ExW — > f/xS'xR be a 

continuous map of the form ijj{x,S^,ri) = (x, 0(x, ^, 77), 77) , with (f){x,(,,ri) positively one- 
homogeneous in (^,r?). The image ^(/x) of any fj, G GY{U;S) is defined as the element of 
GY{U; S') such that 

(/, V'(a^)) (/ ° V-, = (/(a;, ^{x, e, r,), r?), ^(x, ^, r?)) (2.7) 

for every / G G'^°"^ {U xE' xR) . In [6] it is proved that formula (2.7) makes sense also if (f) is 
a Borel map, positively one-homogeneous in (^, r]) , and satisfying \(p{x, ^,r])\ < a|^| +6(a;)|77| 
with a G IR and 6 G L^iU). 

For additional properties of the space GY{U; S) we refer to [6]. 

Systems of generalized Young measures. Let U he a bounded open set in R" , n > 1 , 
and let S be a finite dimensional Hilbert space. If {si, S2, . . . , s„} C {ti,t2, ■ ■ ■ , tm} C M, 
with si < S2 < . . . < s„ and ti < t2 < ■ ■ ■ < tm, the projection 7r*i ;*™ : UxE'^xR —^ 
[/xS"xM is defined by 

Trsi:::',': (.x, 61 , . . . , 6„ , ??) = (a;, 61 , . . . , 6„ , ??) . 

A compatible system of generalized Young measures on U with values in S with time set 
9 C M is a family /j, = {fJ'ti...t^) of generalized Young measures f^ti-.-tm ^ GY{U;E"^) , 
with ti, . . . ,tm running over all finite sequences of elements of with ti < t2 < ■ ■ ■ < tm , 
such that the following compatibility condition holds: 

Msi...s„ = T^l\::*sZ(l^ti...tr^) whenever {si,S2,...,s„} C {ti,i2,...,im}. (2.8) 

The space of all such systems is denoted by SGY{Q, U:E) . 

Given a function p : Q ^ Mb(U; S) , the compatible system of generalized Young measures 
dp associated with p is the element of SGY{Q, U ; S) defined by 

(5p)ti...t„ := ^(p(ti),...,p(t^)) , (2.9) 

where the right-hand side is defined by (2.6). 

Finally, we recall that the variation of /i G SGY{Q, U; E) on a time interval [a, b] , with 
a, 6 G G , is defined as 

k 

Var(/x; a, b) := sup - Htoti...tk {x,^o,---, Cfe, v)) , (2-10) 

i=l 
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where the supremum is taken over all finite families to,ti, . . . ,tk in & such that a = to < 
ti < ■ ■ ■ < tk = b. If /i = (5p for some p: © — > Mb{U ; S) , then Var(jLt; a, b) coincides with 

k 

Var(j>;a,6) := sup^ - , (2.11) 

where the supremum is taken over all finite families to,ti, . . . ,tk in 6 such that a = to < 
ti<---<tk = b. _ 

For the main properties of the space SGY{ld, C/; S) wo refer to [6, Section 7]. 

2.2. Mechanical preliminEiries. We now introduce the mechanical notions used in the 
paper. 

The reference configuration. Throughout the paper the reference configuration O is a 

bounded connected open set in M.'^ , d > 1 , with Lipschitz boundary 90 = Fq U Fi . We 

assume that Fq is closed, H''"^(Fo) > 0, and Fq n Fi = 0. 

On Fo we will prescribe a Dirichlet boundary condition. This will be done by assigning a 

sufficiently regular function w: Fq ^ K.'', or, equivalently, a function w € (Cl-jW^) , whose 

trace on Fq (also denoted by w ) is the prescribed boundary value. 

Every function u G BD{Q) has a trace on dQ, still denoted by u, which belongs to 
(dfl-jW^) . Moreover, there exists a constant C > 0, depending on fl and Fq, such that 

||w||i,n < C \\u\\i,To + C \\Eu\\^^a (2.12) 

(see [30, Proposition 2.4 and Remark 2.5]). 
We will study two alternative situations. 

(DN) Dirichlet- Neumann. In this case we will consider traction- free boundary conditions 
on Fi. 

(DP) Dirichlet-Periodic. In this case O is the cube Q := {—■^, ^)'^ , d>2, and 

To :={xGdQ -.1x11 = 1}, Ti := {x G dQ : \xi\ < 1} , (2.13) 

where x = {xi,x), x = {x2, ■ ■ ■ ,Xd), and we consider only fimctions u which are 
^-periodic, in the sense that u{x + Cj) = u{x) for every x G R"* and every element 
Bi of the canonical basis of M** , with i = 2, . . . ,d. 

Admissible stresses and dissipation. Let K he a closed convex set in M^^^'xR, which 
will play the role of a constraint on the deviatoric part of the stress and on the internal 
variable For every ( £R the set 

K{C) := {a G M^^'^ : (a, C) G K} (2.14) 

is interpreted as the elastic domain and its boundary as the yield surface corresponding 
to ( . We assume that there exist two constants A and B , with < A < B < oo , such that 

{{a, C) G M^^^-^xR : |ap + |Cp < A^} C K C {{a, C) G M^^^'^xM : \a\^ + |Cp < B^} . (2.15) 

We assume in addition that 

{a,C)eK =^ (0,C)ei^. (2.16) 

Let ttr: M^^'^xR ^ M be the projection onto M. The hypotheses on K imply that there 
exist two constants > and 6^ > such that 

MK) = [-aK,bK]- (2.17) 

The support function H: M^^<^xK^ [0,+oo) of K, defined by 

H{^,9):= sup {a:^ + Ce}, (2.18) 
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will play the role of the dissipation density. It turns out that H is convex and positively 
homogeneous of degree one on M^^'^xR. In particular it satisfies the triangle inequality 

if (a + 6, Oi + 62) < H{^i,ei) + H{^2, 62) . 
Prom (2.15) it follows that 

A^\^\^+e^ < H{^, 6) < B./W+^ , (2-19) 
for every {^,6) G M^'^'^xM, and from (2.16) we obtain 

m,.e,>_HM-[;f nil ,2.0) 

for every {^,9) e M^'^'^xM. 

Using the theory of convex functions of measures developed in [13], we introduce the 
functional H: Mb{n;U'j^'^)xMb{Ti) M defined by 



n{p,z) := f_H{^,{x),j^,{x))dX{x), 
Jn 



where A £ M^(il) is any measure such that p << X and z << A (the homogeneity of H 
implies that the integral does not depend on A). Using [13, Theorem 4] and [30, Chapter II, 
Lemma 5.2] we can see that W(p, z) coincides with the integral over fl of the measure studied 
in [30, Chapter II, Section 4], hence H is lower semicontinuous on Mb{Ti;M'j^'^)xMb{Ti) 
with respect to weak* convergence of measures. It follows from the properties of H that Ti 
satisfies the triangle inequality, i.e., 

n{pi+p2,zi + z2)<n{pi,zi) + n{p2,z2) (2.21) 

for every p\,p2& 71^6(0; M^^'^) and every Zi,Z2e MbiP). 

For every e > we introduce the function : M^^'^xR R defined as 

H,{^,e) :=7J(^,e) + ||e|2 + ||(?|2^ (2.22) 
and the corresponding integral functional Hs: L'^{Vt;M'^j^'^)xL'^{^) M defined by 

T-(-£{p,z)~ / Hsip{x),z{x))dx. 
Jn 

The convex conjugate iJ* : M'^'^^'xR -^R of is defined by 

H:{a,0~ sup {a:^ + C9-H,{^,e)}. 

(C,f)6M^^'*xK 

Since the convex conjugate H* of H satisfies H*{a,() = for (tr, C) € K and H*{a,Q = 
+0O for (ct, C) ^ K (see [27, Theorem 13.2]), using [27, Theorem 16.4] one can prove that 

H:iaX) = i-JiaX)~PK(a,C)\\ (2.23) 

where Pk ■ M^^'^xR ^ ii' is the projection onto K. This implies that H* is differentiable, 
and that its gradient is given by 

N^a, C) := I {{a, Q - Pxia, Q) • (2.24) 

Note that Nj^ is Lipschitz continuous. 

Let H* : L^(f2; M^^'') x L'^{n) -» R be the convex conjugate of He - By a general property 
of integral functionals (see, e.g., [10, Proposition IX. 2.1]) we have 



KiaX)^ I H;i<j{x),C{x))dx, 
Jn 



so that, by the Dominated Convergence Theorem, its gradient is given by 

dn*{a,C)ix) = Nl({a{x),Cix)), for a.e. x G 17 . (2.25) 
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Therefore dH.* is Lipschitz continuous. 

The elasticity tensor. Let C be the elasticity tensor, considered as a symmetric positive 
definite hnear operator C: Mf^^f^ ^iym ■ We assume that the orthogonal subspaces M^^** 
and M/ are invariant under C. This is equivalent to saying that there exist a symmetric 
positive definite linear operator and a constant k > such that 

C^:=Ci5^D + K(trO/ (2.26) 

for every ^ e M^^^ . Note that when C is isotropic, we have 

Ce = 2/x^z? +K(tre)J, (2.27) 

where /z > is the shear modulus and k > is the modulus of compression, so that our 
assumptions are satisfied. 

Let Q: M^^^ [0, +oo) bo the quadratic form associated with C, defined by 

0(0 := ice : e = \^D^D ■■ + f (trO' • (2.28) 
It turns out that there exist two constants etc and /3c , with Q < ac < 13c < +oo , such that 

aclCP < QiO < Pc\e (2.29) 
for every ^ G ^tym ■ These inequalities imply 

|Ce| < 2(3c\^\ . (2.30) 

The softening potentiaL Let F : M ^ M be a function of class , which will control the 
evolution of the internal variable Cj and consequently of the set K{() of admissible stresses. 
We assume that there exists a constant M > such that 

-M < V"{e) < for every 6 , (2.31) 
-bK < V'{+oo) < V'{-oo) < UK , (2.32) 

where V'{zkoo) denote the limits of V'{9) as 6 ^ ±00, while qk and are the constants 
in (2.17). We denote the recession function of V by V°° ; it is defined as 

V-ie):= lim M = ^;(-oo)^ if^^O^ (2.33) 
^' t [V'(+oo)6' if6i>0. ^ ' 

Let us fix £ e (0, 1) such that 

V {-00) < (1- e)aK and V {+00) > -{1 - e)bK ■ (2.34) 

Then, by (2.20), 

H{^2-^i,e2-ei) + v{e2)-v{ei) > 

> eH{^2 - ?i, ^2 - ^1) + (1 - e)H{Q,e2 - 9{) + V'{6i){62 - ^1) 

for a suitable 6*3. By (2.20) and (2.34), using the monotonicity of V , we obtain (1 — 

e)H{id,e2-ei) + V'{ez){e2-ei) >0, hence 

H{^2 - 6, 02 - Oi) + V{92) - Vi0i) > eH{^2 - €1, ^2 - ^1) . 
Therefore (2.19) implies that there exists a constant Cy > such that 

^^(6 -^1,02- Oi) + ¥(62) - V{9i) > 16 -ii\ + C^ \02 - 01 1 (2.35) 

for every ^1,^2 & M'j^'' and every 6*1, 6*2 € M. 

It is convenient to introduce the function V : (f2) K defined by 

V(^) := / V{z{x)) dx 
Jq 

for every z £ L-^ (Q) . 
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The prescribed boundary displacements. For every t G [0, +00) we prescribe a bound- 
ary displacement w{t) in the space H^{Cl; M'*) . This choice is motivated by the fact that we 
do not want to impose "discontinuous" boundary data, so that, if the displacement develops 
sharp discontinuities, this is due to energy minimization. 

We assume also that w G ^/^^([O, +00); il^(f2; K"^)) , which means, by definition, that 
w G H'^{[0,T]; H^{Q;R'^)) for every T > 0, so that the time derivative w belongs to 
L^{[0,T];H\n;W^)) and its strain Ew belongs to L^{[0,T]; L^iniM^^^)) . In the Dirichlet- 
Periodic case (DP) we assume also that w{t) is (the restriction to Q of) an ^-periodic 
function. For the main properties of Sobolev functions with values in reflexive Banach 
spaces we refer to [2, Appendix]. 

Elastic and plastic strains. Given a displacement u G BD{^l) and a boundary datum 

w e H'^in; M'') , the elastic strain e e L'^{n; Mf^^) and the plastic strain p e Mb{Tt; M^""^) 
satisfy the weak kinematic admissibility conditions 

Eu = e+p inf7, (2.36) 

p= {w-u)Qnn'^-^ onFo, (2.37) 

where n denotes the outward unit normal. The condition on Fq shows, in particular, that 
the prescribed boundary condition w is not attained on Fq whenever a plastic slip occurs at 
the boundary. It follows from (2.36) and (2.37) that e = E°-u — p^ a.e. in O and = E^u 
in O. Since ivp = 0, it follows from (2.36) that divw = tre e L'^{Q) and from (2.37) that 
{w-u)-n = Q H'^-i-a.e. on Fq. This shows that, if u e H'^ {Q.;m.<^) and p e L'^{Q.;M'^j}'^) , 
a stronger kinematic admissibility condition is satisfied. 

Definition 2.1. Given w € H^{^;W^), the set A{w) of admissible displacements and 
strains for the boundary datum w on Fq is defined as the set of all triples (w, e,p) , with u € 
ili(0;R''), e e L2(0;M^^4), p e L'^ {Q.^Mj}''-) , which satisfy the kinematic admissibility 
conditions 

Eu = e+p a.e. in r^, (2.38) 
u = w 7i:'^-i-a.e. on Fq . (2.39) 

In the Dirichlet-Periodic case (DP) we always assume that w is the restriction to Q of an 
^-periodic function belonging to Hlg^{M.'^;M.'^) and we add the requirement that the same 
property holds for u . 

The stress. The stress a G 1.2(0; M^^^) is given by 

C7 := Ce = Ccez) + Ktre, (2.40) 

and the stored elastic energy by 

Q(e)= / Q{e{x))dx = ^{a,e). (2.41) 
Jn 



It is well known that Q is lower semicontinuous on L'^{Cl; M^y^) with respect to weak 



convergence 

sym 



Let Q* : L'^{Q;M'^y^) [0, +00) bo the convex conjugate of Q. It is well known that 



Q*(a) = / Q*{a{x))dx, 
Jn 



where Q* : Mf^^ ^ M is the convex conjugate of Q given by 

Q*(<t) = ia:C-V. (2.42) 

If cr G 1/2(0; Mf^^) and diver e L'^{Q;R'^), then the trace of the normal component of 
a on dfl , denoted by [an] , is defined as the distribution on dfl such that 

(M,^>an := (diva, V) + {(T,Ei;) (2.43) 
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for every S H^{n;R'^). It turns out that [an] G H-'^^^{dn;R'^) (see, e.g., [30, Theo- 
rem 1.2, Chapter I]). We say that [an] = on Ti if {[an],tp)dn = for every ip e H^{n;R'^) 
with ip = W^-^-a.e. on Tq. 

3. The minimum problem 

In this section we study in detail the incremental minimum problems used in the discrete- 
time formulation of the regularized evolution. The data arc the current values pq <E 
L^[fl;M.'^'^) and Zq G L-^(ri) of the plastic strain and of the internal variable and the 
updated value wi G H^{fl;R'^) of the boundary displacement. By solving the minimum 
problem 

min{Q(e) + n{p -po,z- zq) + V(^) + - Po||i + zo]\l} (3.1) 

over the set of all {u,e,p) E A{wi) and all z G L'^(n), we get the updated values w, e, p, 
and z of displacement, elastic strain, plastic strain, and internal variable. In (3.1) £ > is 
a prescribed "viscosity" parameter, while r > will play the role of the time discretization 
step. 

We are now in a position to prove the existence of a solution to (3.1), provided r < s/M , 
where M is the constant appearing in (2.31). 

Theorem 3.1. Let e > 0, t £ {0,e/M), wi G H^{n;R'^), po G L^{n;M'j^''') , and zq G 
L^{Cl). Then the minimum problem (3.1) has a unique solution. 

Proof. Let {uk,ek,Pk,Zk) be a minimizing sequence. By (2.35) we have 

'H{pk -po,Zk - zo) + V{zk) > Cy]\pk -po\\i + Cy\\zk - 2;o||i + V(2:o) ; 

therefore, (3.1) implies that Q{ek) + ^Ibfe —Polli + ^Ikfc ^ ^olli is bounded. Consequently, 
the sequences pk and Zk are bounded in L^{n; M'^^'') and in L^{Q) , respectively, and Q{ek) 
is bounded. Using (2.29) we deduce that the sequence is bounded in L^(r2; Mf^^) . 
Since Euk = Ck + Pk a-e. in fl, it follows that Euk is bounded in i^(f2;Mfy^). Since 
Uk = wi TC'^~^-a..c on Fq, Uk is bounded in H^{n;R'^) by the Korn-Poincarc inequality 
(see e.g. [4, Theorem 6.3-4]). Since r < e/M , the functional in (3.1) is convex, hence weakly 
lower semicontinuous. The existence of a minimizer follows now from the direct methods 
of the calculus of variations. The uniqueness is a consequence of the strict convexity of the 
functional. □ 

We now derive the Euler conditions for a minimizer of (3.1). In what follows, the symbol 
d denotes the subdifferential of convex analysis ( sec, e.g., [10, Definition 1.5.1]). 

Theorem 3.2. Lei e > 0, r G (0,£/M), po G L^{n;M'^j^'^) , Co S L^{^), wi G H\Q;R'^) , 
(«i,ei,pi) G A{wi), zi G L\n), ai := Cei G ^^(OjM^;^), and Ci := -V'{z^) G i°°(0). 
// {u\,ei,pi,z\) is the solution of (3.1), then 

W(P +Pl-P0,z + Zi- Zq)- nipi -po,Zi- Zo) > 

> -(^1,?:) + (Ci,5) - f (pi -po,p) - f (zi - zo,z) 

for every {u,e,p) G ^(0) and every z G L'^i fl). Condition (3.2) implies that 

divo"! =0 m il , [ain] =0 onTi , (3-3) 

((o-i)r> - 7(pi -l'o),Ci - f (^1 - ^o)) e dH{pi -_po, zi - zo) . (3.4) 

In the Dirichlet- Periodic case (DP) the condition [an] =0 on Fi is replaced by the re- 
quirement that a is the restriction to Q of an x -periodic function satisfying diver = 

in R"^. 

Conversely, conditions (3.3) and (3.4) imply that (ui, ei,pi, Zi) is a solution of (3.1). 
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Proof. To prove (3.2) we fix {u,e,p) S A{0) and z G L^(f2). For every s > the triple 
(ui + su, ei + se,pi + sp) belongs to A{w-i) , and hence by minimality, 

Q{ei + se) ^-^{{sp + pi - po,sz + zi - zo) + V{z\ + sz) + 

+^||sp + Pi -Polli + ^\\sz + zi - zoWl > 

> Q(ei) + Hipi - po, zi - zo) + V{zi) + ^Ibi - PoWl + - zo\\l . 

Using the convexity of H we obtain 

s{U{p + pi -po,z + Zi- zo) - H{pi -po, zi - zo)} > 
> -Q(ei + se) + Q(ei) - V(zi + sz) + V(zi) - 
- ^Ibi + sp-po\\l + ^Ibi -Polli - if\\zi +SZ- zoWl + ^\\zi - zoWl 

for every s > 0. Taking the derivative with respect to s at s = 0, we obtain (3.2). 

Taking e = Eu, p = 0, and z = in (3.2) we deduce that {ai,Eu) = for every 
u G H^{Q) with u = H'^~^-a..e on Fq, which gives (3.3) thanks to (2.43). The changes in 
the Dirichlet-Periodic case are obvious. 

Taking '0 = and e = —p, we obtain from (3.2) 

Hip + pi-po,z + zi- Zo) - H{pi - Po, zi - Zo) > 

> {cFi.p) + (Ci,z) - 7(Pi -POiP) - 7(21 - zo,z) = 
= {{<^i)d,p) + {Ci,z) - ^{pi -PQ,p) - 7(21 - za,z) 

for every p G L'^{n;M'j^'^) and every z e L'^i^l), which implies (3.4). 

The fact that (3.3) and (3.4) imply minimality follows from the strict convexity of the 
functional. □ 

The following theorem provides a dual formulation of problem (3.1). 

Theorem 3.3. Let e > 0, t € {0,e/M), wo € H^{Q;R''), {uo,eo,po) € A{wo) , zo G 
L^{fl), (To := Ceo, O'^d Co •= —y'izo)- Let Wi G iJ-^(f2; R'') , let {ui,ei,pi, Zi) be the 
solution of (3.1), and let ai := Cei and Ci := —V'{zi). Then 

{^{pi-po),H^i-zo))=dniiiai)D,Ci)- (3.5) 

Moreover {ai,(i) is a solution to 

min{iQ*(a - ao)+ntiaD,0 - i(C,^i - ^0) - ^{(t,Ewi - Ewo)} (3.6) 

over the set of all {a,Q & -^^(f^;Mf^^)xZ/^(n) with diva = in fl and [an] = on 

Fi . In the Dirichlet-Periodic case (DP) the condition [an] =0 on Fi is replaced by the 
requirement that a is the restriction to Q of an x -periodic function satisfying diva = 
in M."^. 

Proof. As dH is positively homogeneous of degree 0, formula (3.4) can be written as 

{{(Ti)d - f (Pi -Po), Ci - H^i - zo)) e aW(i(pi -po), i(zi - zo)) , 
which is equivalent to 

((c7l)z3,Cl) e dHei^ipi - Po), H^l - zo)) . 

By a general duality formula (see, e.g., [10, Corollary 1.5.2]) this is equivalent to (3.5). 
Using (2.42), we obtain 

7Q*(fT-ao)> iQ*(ai-ao) + i(a-ai,ei-eo) (3.7) 
for every a € L2(f]; Mf^^„^J . By (3.5) we have 

nt{aD,0 >n:{{ai)D,Ci) + H<^D-{(Ji)D,Pi~Po) + HC-Ci,zi- zo) = 
= K{{(^i)d, Ci) + 7(^ - ^I'Pi - Po) + 7(C - Ci, ^1 - ^0) . 
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(3.8) 



(3.9) 



Since po = Euq — cq and pi ~ Eui — ei , wc obtain 

H*(aD, C) >ni{{ai)D,(:i) + \{<J-<JuEui- Euq) - 
- 7(0-- cri,ei - eo) + ^(C- Ci,^i - 2^0) 

for every ((7,0 e L'^{Q.;M'l^^)xL'^{Q) . Adding (3.7) and (3.8) we deduce 

^Q*{<y- To) + UHudX) - 7(C, ^1 - -zo) - H'^.Eui - Euo) > 
> iQ*(ai - ao) +H*((ai)D,Ci) - 7(Ci,^i - ^o) - ^{^uEui - Euo) . 

If div (T = in f2 and [an] = on Fi , integrating by parts and using the equality Ui — uq = 
wi — Wo W^~"^-a.e. on Fq, we obtain {Eui — Euo,cr) = {Ewi — Ewo^a). Together with 
(3.9), this proves that (ci, Ci) is a solution to (3.6). □ 

4. Regularized evolution 
In this section we study the notion of regularized evolution. 
4.1. Definition and properties. We begin with the definition. 

Definition 4.1. Let w e Hl^{[{),+oo)]H^{n;M.^)) , Uq G H^in;R'^), eo € ^^(f]; Mf^^^) , 
Po € L^(ri; M^^'') , zq e L^(ri), and let £ > 0. A solution of the e -regularized evolu- 
tion problem, with boundary datum w and initial condition {uo,eo,Po, Zo) is a function 

iUe,ee,P^,Ze), with 



We e ifi„,([0,+oo);iJi(12;M^)), G 7?/„,([0, +00); ^^(r!; Mf-^)) , 

Pe e /fL([0, +00); L^(f2; M^x^^)) , e if/„,([0, +00); L\n)) , 



(4.1) 

such that, setting 

the following conditions arc satisfied: 

(evO)^ initial condition: (itg (0) , (0) , (0) , Zg (0) ) = (wq , eo , po , 2:0 ) ; 
(evl)^ kinematic admissibility: for every t G (0, +00) 

Eus{t) = ee{t) -\- p^{t) a.e. inO, 
Us{t) = w{t) n^-^-a.e. in Fq ; 

(ev2)£ equilibrium condition: for every t G (0, +00) 

divcTsit) = in O , [a-e{t)n\ = on Fi ; 

(ev3)e regularized flow rule: for a.e. t G (0, +00) 

{p^{t),z,{t)) = iV^(a,(0z3,CW) a.e. in O, 

where JV|. is defined by (2.24). 

In the Dirichlct-Periodic case (DP) we assume that for every t G [0, +00) the functions 
w{t), Mo, and u^{t) are restrictions to Q of S-periodic functions of class i?;Q^(R''; R'') , 
and in the equilibrium condition (cv2)g the equality [cr£(t)n] = on Fi is replaced by the 
requirement that erg (t) is the restriction of an x-periodic function with div eXe (t) = in M'* . 

Remark 4.2. Let us fix t > such that the derivatives Pg{t) and Zs{t) exist. Then the 
following conditions are equivalent: 

{p,{t), ze(t)) = Nf,{crS)D, CM) a.e. in ft , (4.2) 
{crS)D, Ceit)) e m^iPM, ^S)) a.e. in O , (4.3) 
{<T,{t)D - ePe{t), - ^zS)) G dn{p,{t), z,{t)) a.e. in O . (4.4) 

Indeed, by (2.25) wc have 5W*((Te(t), C|e(t)) = N 1^ {a ^ {t), Ce {'<'')) ^ so that (4.2) and (4.3) are 
equivalent by a standard property of conjugate functions (see, e.g., [10, Corollary 1.5.2]). 
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The equivalence between (4.3) and (4.4) follows immediately from the definition of He 
(see (2.22)). 

The following theorem shows that the modified fiow rule (ev3)e can be replaced by a suit- 
able stress constraint and an energy equality. To formulate these conditions it is convenient 

to introduce the convex set 

K-ifl) -.^ {ia,C) e L^{n;m'j^'^)xL°°{n) : {a{x)X{x)) e K for a.c. .t G 0} . (4.5) 

Theorem 4.3. Let w , uq, cq , pq , zq, and e be as in Definition 4-1, and let (it^, e^, p^, 2;^) 
be a function satisfying (4.1). Then {ue,ee,p^, z,.) is a solution of the s -regularized evo- 
lution problem with boundary datum w and initial condition {uq, eo,Po, zq) if and only if 
it satisfies the initial condition (evO)e, the kinematical admissibility (evl)e, the equilibrium 
condition ( ev2)e , and the following properties: 

(ev3)e modified stress constraint: for a.e. t G (0, +oo) 

(eTe(i)B - spM, cm - ^zS)) e /C(fi) ; 



(ev4)£ energy equality: for every T > we have 

Q{ee{T))+ r n{p,{t),zS))dt + V{ze{T)) + e C \\p,{t)\\ldt + e C \\zS)\\ldt = 
Jo Jo Jo 

= Q(eo) + V(zo) + / {a,{t),Ew{t))dt. 
Jo 

Proof. Suppose that {ue,ee,Pe, z^) satisfies (evl)e, (ev2)e, and (ev3)e. As H is positively 
homogeneous of degree one and dH{0, 0) = K (see, e.g., [27, Corollary 23.5.3]), by a general 
property of integral functional (see, e.g., [10, Proposition IX.2.1]), we have 

dn{p'{t),z'{t)) c dn{o,o) = icin) . 

Therefore (4.4) implies (ev3)e . 

Since H is positively homogeneous of degree one, the Euler relation gives {a, p) + (C, z) = 
Ti,{p,z) whenever (u, C) G d'H{p,z). Therefore, (4.4) implies 

n{p,it), Ze{t)) = {aeit)D - ep,{t),p,it)) + (Ce(t) - eze{t),ze{t)) , (4.6) 
which is equivalent to 

HiPM, ^e{t)) + {V'{zS)),Ze{t)) + e\\pM\l + eW^eiml = {cTe{t),pM) ■ (4-7) 
By (evl)£ we have 

{ae{t),p,{t)) = {ae{t),Eu,{t)) - {ae{t),ee{t)) . (4.8) 
Integrating by parts and using (ev2)£ , we obtain 

(<T,(t), EiiS)) = {(Teit), Ew{t)) . (4.9) 
Combining (4.7), (4.8), and (4.9), we deduce that 

{aeit),ee{t))+n{p,{t),ze{t)) + {V'{ze{t)),ze{t)) + 

+ e\\pM\l+4^s{tm-{'re{t),Ew{t)). ^ ■ ' 

The energy equality (ev4)e can be obtained from (4.10) by integration. 

Conversely, assume that {us,es,p^,Ze) satisfies conditions (evl)^, (ev2)£, (ev3)e, and 
(cv4)£. By differentiating (cv4)e we obtain (4.10). Thanks to (4.8) and (4.9), from (4.10) 
we deduce (4.7), which is equivalent to (4.6). By (evS)^ for a.c. t G (0,+oo) we have 

ia,{t)D - ep.it), CM - ez,{t)) e - 5W(0, 0) . (4.11) 
Using the homogeneity of ?i, condition (4.4) follows easily from (4.6) and (4.11). □ 
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Theorem 4.4. Let w, uq, cq, Po, zq, and e be as in Definition ^.1, and let ctq := Ceo- 
Suppose that the following conditions are satisfied: 

(evl)o kinematic admissibility: 

Euo = cq+Po a-e. in fl , 
uo =w{0) H'^-'^-a.e. m Tq . 



(ev2)o equilibrium condition: 

div ao = in fi , [uon] ~ Q onVi . 

Then there exists a unique solution of the e -regularized evolution problem with boundary 
datum w and initial condition {uo,eo,po,zo) ■ 

The proof is postponed to the next subsection. 

Let M''^'' be the space of all dxrf-rnatriccs. The uniqueness result of the previous theo- 
rem allows to simplify the problem in the spatially homogeneous case, when the boundary 
condition has the form 

w{t, x) := ^{t)x , To := 90 , (4.12) 

with 

?effi,([0,+oo);M'^x''), (4.13) 
while the initial condition has the form 

uo{x) := ^ox , eo(x) := , po{x) := , zo{x) := 6*0 , (4.14) 

with 

^0 e M'*^'' , ^0 e M^ym , Co e M^^'' , 6»o e M . (4.15) 
We assume the following compatibility condition 

Co -eo+a, m=io, (4.16) 

where denotes the symmetric part of • 

Proposition 4.5. Assume that w , Tq , uq , cq , Po , zq satisfy (4.12)-(4.16), and let e > 0. 
Then {us,ee,Pg, z^) is the solution of the e -regularized evolution problem with boundary 
datum w and initial condition {uo,eo,po,zo) if and only if 

uS,x)=mx, ee{t,x)=^l{t), pS,x)=eAt), zS,x)=e,{t), (4.17) 

with 

m--=^%t)-es{t), (4.18) 

where ^^{t) denotes the symmetric part of ^{t) and (C|,^e) is the unique solution in 
^^(oc([0,+oo);Mf^>=^)x//i„^([0,+(X3)) of the Cauchy problem 

{t{t)-m,^e{t)) = Ni,{CnmD,-v'{e,m, (4.19) 

^1(0) = eM = eo. (4.20) 

Proof. Let {(,^,6^) be the solution of the Cauchy problem (4.19), (4.20), whose existence 
follows from the Lipschitz-continuity of and V . By the definition of Nf^ (see (2.24)) 
it follows that t{t) - £,l{t) € M^'''' for a.e. t e [0,+oo). Moreover, by (4.15), (4.16), 
and (4.20) we have also ^'(0) - ^1(0) e M^'^''. It follows that ^(i) - ^tit) G M^'^'' for 
every t e [0, +oo), hence ^|?(t) e M^^'^ by (4.18). It is then easy to see that the function 
(tie, ee,Pg, Zg) defined by (4.17) satisfies all conditions of Definition 4.1. □ 
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4.2. Proof of the existence theorem. To prove Theorem 4.4 it is useful to introduce 
the notion of iJ-dissipation of a possibly discontinuous pair of functions p: [0,+oo) — > 
Mb{Tl;M'j^'') and z: [0, +oo) Mb{Ti) . For every pair of times a, b € [0, +oo), with 
a <b, the ff -dissipation of {p,z) on the interval [a,b] is defined by 

k 

VHip,z;a,b) :=sup5^W(p(t,) -pfe_i),z(i,) -zfe_i)), (4.21) 

i=i 

where the supremum is taken over all finite families to,ti, . . . ,tk with a = io < ^i < • • • < 
tk = b. For every p e /f/„';!([0, +oo); M^^"^)) and every z e hII{[0,+oo); L^{n)) we 

have ^ 

VH{p,z;a,b)= f H{p{t),z{t))dt (4.22) 

J a 

(see, e.g., [5, Theorem 7.1]). 

Proof of Theorem 4-4- Wc prove the existence of a solution to the ^-regularized evolution 
problem by time discretization, using an implicit Euler scheme, which leads to an incremental 
minimization problem. The use of the dual formulation of the problem is inspired by [29]. 
The proof is divided into several steps. 

Step 1. The incremental problems. Let us fix a sequence of subdivisions {t].)i>Q of the 
half-line [0,-|-cxd), with 

Q = tl<ti<---<f-^ <tl^+<^, (4.23) 
Tfe := sup(4 - ^ asfc^oo. (4.24) 

i 

It is not restrictive to assume also that < e/M for every k, where M is the constant 
introduced in (2.31). 

For every i we set u>^ := w(t].) . Wc define e\, and by induction on i. We 
set (Wfc, e^,p^, 2:^) := (uq, eo,po, -^o) , and for i > 1 we define {u\,e\,p\, z^) as the solution 
to the incremental problem 

min{Q(e) + n{p - p^-^z - z^-') + V(^) + ^A\p- pf'lli + 2%ll^ " 4"'ll2} . (4-25) 

where := t], — t]^^ and the minimum is taken under the conditions {u,e,p) G ^(w^) 
and z e L^(f]). The existence and uniqueness of the solution to this problem is proved in 
Theorem 3.1. 

For every i > we set CTj. := Ce^. and Q := —V'{zl). We consider the piecewise constant 
interpolations defined by 



Mfe(t) := ui , ek{t) ~ el , p^{t) , Zk{t) := zi , 

<^k(t):=cjl, Ckit)-=a, Wk{t):=wl, [t]k:=tl 



(4.26) 



for t G [tl,tl+^). By definition 

{uk{t),ek{t),Pk{t)) e A{wk{t)) (4.27) 

and 

Uk{0)=UQ, efc(0) = eo, Pfe(0)=jJo, Zk{0) = ZQ. (4.28) 

A solution of the e -regularized evolution problem will be obtained by taking the limit, as 
fc — > oo, of these piecewise constant interpolations. We will also consider the piecewise affine 
interpolations 



Pkit) 



= el + {t-tl){el+'^ei)/Tl+\ 

= Pi + it-ti)ipl+'-pl)/rl+\ (4.29) 



4 + (i-4)(^^'-4)K+S 
= 4 + (i-4)(^r'-4)/4+' 
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for every tG[tltl+'). 

We derive now an energy estimate for the solutions of the incremental problems. 

Lemma 4.6. For every T > there exists a sequence — > 0+ such that 

Q{ek{t2)) + Vh{pu, Zk-^tiM) + V{zk{t2)) + 



+ 1 r^''\\i>t{mldt+'^ r^^''\\ztmldt< (4.30) 



Q{ek{ti))+V{zk{ti))+ t ' {ak{t),Ew{t))dt + iol . 



'[ti]fc ^ J[t 

\t2] 

< Q{ek{ti))+V{zk{ti)) + 

'[ti\k 

for every k and every t\,t2 € [0, T] with t\ <t2- 

Proof. Let us fix T > 0. Since (p^., Zk) is piecewise constant on [0, T] , it is easy to see that 
VH{Pk,Zk;t„t2)= nPl-Pr\4-4-'), (4.31) 

r=i+l 

where i and j are such that t]. = [ti]k and tj. = [t2]k- Therefore, we have to prove that 
there exists a sequence 0+ such that 

Qi4)+ E n{pi-pr\zi-zr')+vi4) + 

r—i-'rl 

+ I r \\ptmldt+l f" Wztmldt < (4.32) 
< Qiel) + V{zl) + / {ak{t), Ew{t)) dt + u^. 

for every i, j, k with < z < j and t\<T . 

Let us fix an integer r with i < r < j and let u := u^" — w^r + w^. and e := 

el~^ - Ewl~^ + Ewl; then, {u,e,plr^) G Testing with (u, e,p^"\ 4"^) , by the 
minimality condition (4.25) we have 

Q{el) + nipl - pr\ zl - zl-') + V{zl) + 

+ 2%lbfe - pr'lli + 2%I14 - < (4.33) 

<Q(er^+£;i«^-£;i«r^)+v(zr'), 

where the quadratic form in the right-hand side can bo developed as 

Qie^r' + Ewl - Ewl-') = Qiel-') + {al-\Ewl Ewlr') + Q{Ewl - Ew^-') . (4.34) 
From the absolute continuity of w with respect to t wc obtain 



^ .)dt, 

"k 

where the right-hand side is a Bochner integral of a function with values in H'{Q.; W^) . This 
implies that 

Ewl-Ewl-''=l Ew{t)dt, (4.35) 

where the right-hand side is a Bochner integral of a function with values in L^(0;M^^^). 
By (2.29) and (4.35) we get 

Q{Ewl - Ewl-') < \\Ew{t)hdty. (4.36) 
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(4.37) 



Note that we can also write 

ll4-4-'ll2 = (*;-C') f\\\Kmldt. 

By (4.33)-(4.37) we obtain 

Q{ei) + n{pi - p^^ 4 - ^r') + v(4) + 1 fj_^ m mi + ii^^ mi) dt < 

< Q{el-')+V{zr')+ f\{al-\Ew{t)}dt + Pc( f\ \\Ew{t)hdtf < (4.38) 
<Q{er')+V{zr')+ f\{ar\Ew{t))dt + pl f\\\Ew(t)hdt, 



where 

i-tl 



by the absolute continuity of the integral. Iterating now inequality (4.38) for i + 1 <r < j , 



we get (4.32) with lj^ := pi \\Ew{t)\\2 dt. □ 

We now prove a dual energy estimate, where ac and M are the constants introduced in 
(2.29) and (2.31). 



Lemma 4.7. We have 

r-[t: 
'[tih 

-[t2]k fit 

< 

'[*l]fc " -^Itl 



ac [ ^ \\e^mldt + K{ak{t2)D,Ckit2)) -K{<7k{h)n,Ckiti)) < 

M \\ztm\ldt+ t'^^\&t{t),Ew{t))dt, 

■J[ti]k •>[ti]k 



(4.39) 



for every k and every ti,t2 € [0, +oo) with ti <t2- 
Proof. We have to prove that 



«C t^\\etmldt + n:{{ai)n,a)-K{{al)n,a))< 

<M [ " \\z^{t)\\ldt+ [ \&t{t),Ew{t))dt, 
Jti Jti 



(4.40) 



for every k and every < i < j . Let us fix an integer r with i < r < j . Wc observe that 
divcr^^^ = in ri and [cr^~^?^] = on Fi . This follows from (3.3) if r > 1, and from the 
equilibrium condition (ev2)o if r = 1 . By Theorem 3.3 we have 

^Q*K - <^r') + ni{{ai)nXl)~ni{{ai-^)u,cr') < 

< - Cu-\zl - zl-') + Mai - ar\Ewl - Ew^^') . 
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Since ac\\el-el~'\\l < Q*{al-a'^') by (2.29) and (C,--Cr\4-4"') < 
by (2.31), we obtain 

"c r \\Kmidt+n;{{ai)n,a)-K{{ai-')D,a-')< 

<M f\ \\z^mldt+ r {cTt{t),Ew{t))dt. 



k 



Summing over r we obtain (4.40). □ 

Proof of Theorem 4-4- Step 2. Proof of the bounds. Let us prove now that for every T > 
there exists a constant Ct , independent of k and e , such that 

sup \\ek{t)h < Ct , Vnip^, Zk;0,T) < Ct , (4.41) 
te[o,T] 

f \\Kmldt<-^, f \\p-{t)\\ldt<^, f \\zt{t)\\ldt<^. (4.42) 
Jo t '\ i Jo £ Jo ^ 

By (2.35) and (4.21) for every t € [0, +oo) we have 

VuiVk^Zk-M + V{zk{t)) >n{Pk{t) - Po.Zk{t) - zo) + V{zk{t)) > 

>C^\\p^{t)-po\\i + C^\\zk{t)-zoh+V{z,). > 



Let us fix T > 0. Using the discrete energy inequahty (4.30) with = and t2 = t < T 
and inequahties (2.29), (2.30), and (4.43), we deduce that 



ac\\ekml<Pc\\eo\\l + 2Pc sup \\ek{t)\\2 [ \\Ew{t)\\2dt 

te[o,T] Jo 



for every k and every t G [0, T] . The first estimate in (4.41) can be obtained now by using 

the Cauchy inequahty. 

By (4.30) and the first inequahty in (4.41) we have that 

VHiPk,Zk;0,t)+V{zk{t))<C (4.44) 

with C independent of k and t. By (4.43) this implies the boundedness of ||P/c(i)l|i ^'^d 
||zfc(t)||i, and in turn, the boundedness of V{zk{t)). Now we can use (4.44) to obtain the 
second estimate in (4.41). The last two inequahties in (4.42) foUow immediately from (4.30) 
and (4.43). The first inequality in (4.42) can be obtained from (4.39) with ti = and 
t2=T thanks to (2.23) and (4.28). 

To continue the proof of Theorem 4.4 we need the following lemma, based on Gronwall's 
inequality. 

Lemma 4.8. For every T > the sequences Cfe and Zk satisfy the Cauchy condition in 
L°°([0,T];L2(r2;Mfy^^)) and L°°{[0,T]; L'^{Q)) , respectively. 

Proof. By (2.25) and (3.5) for every i and k we have 

{plzl)-{pl-\zr')=TlN^^{{al)o,-V'{zl)). 
Since Eul = e\+p\ and Eu'iT^ = e^~"^ + p'l^^ ^i-C- in Q, wc obtain 

i-el zl) i-el-\ zl7') = TlN^A{<^l)n, -V'{zl)) - {Ev^^, 0) + {Ev^^-\Q) . 
Summing for 1 < i < j we obtain 

(-4,4) - (-eo,^o) = Y.riNl,{{ai)n, -V'{zi)) - {Eulo) + {Euo,0) . 
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Since TV^ is 1/e-Lipschitz and V is M-Lipschitz, we have 
(-4,4) - i-eo,zo) = i2^iNk{{ai-')n,-V'{zi-')) + Ri - (^4,0) + {Euo,0) , (4.45) 

i=l 

where the rest Rl G L^{Q;M'j<^'^)xL'^{n) can be estimated by 

iii?iii2 < li2^i{\wi - 4-% + M 114 - zi-%) < 



<^ [\2Pc\\etm2 + M\\ztm2)dt, 
£ Jo 



thanks to (2.30). Using (4.42) we can prove that for every T > there exists a constant 
A"^ , depending on e , but independent of j and k , such that 

\\Rih<A'^Tk (4.46) 

for every j and k with tj, <T . 

If t € [0, T] , for every k there exists a unique j such that t-j, <t < t-j^^ . From (4.45) we 
obtain ^ 

(-efc(t), Zk{t)) - i-eo, zo) = [ NI^{<Tk{s)D, -V'{zk{s))) ds + 

•^0 (4.47) 
+ Rk{t) - {Euk{t),0) + {Euo,0) , 

where ^ 

Rk{t)=Ri- I NI,{ak{s)D,-V'{zk{s)))ds. 

In these formulas we are using Bochncr integrals of functions with values in i^(f2;M^^'^) x 
1,2(0). It follows from (2.32), (2.24), (4.41), and (4.46) that for every T > there exists a 
constant , depending on e , but independent of t and k , such that 

\\Rkit)h < B'^Tk (4.48) 

for every t € [0, T] and every k. 

Wc now consider two indices h and k. Subtracting term by term the equations corre- 
sponding to (4.47) wc obtain 

(-efc(t), ^fc(t)) - {-eh{t),Zh{t)) = Ihk{t) + Rk{t) - Rh{t) + 
+ {Ewk{t) - Euk{t),0) - (Ewhit) - Euhit), 0) - {Ewk{t) - Ewh{t),0) , ^^"^^^ 

where 

Ihkit) := f {NI,{ak{s)D, -V'{zk{s))) - NI,{<ju{s)d, -V'izhis)))} ds . 

Since Nj^ is 1/e-Lipschitz and V is M-Lipschitz, using (2.30) we obtain that there exists 
a constant L, depending on e, but independent of t, /i, fc, and T, such that 

\\Ihk{t)\\2 <L f iphk{s)ds (4.50) 
Jo 

where 

^hk{t) := {\\ek{t) - Bhiml + M) - ^ft(t)lli}'^' • (4.51) 

Since w G H^{[0,T]; H^{rt;M.'^)) , there exists a constant , independent of t and k, such 
that \\Ew{t) - Ew{s)\\2 < C^jt- sj^/^ for every s,t e [0,r]. It follows that 

WEwkit) - Ew,{t)h < C^ir'J'+r'J') (4.52) 
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for every t € [0,T] and for every h and k. 

Since div{cTk{t) — o-h{i)) = in O and [(o'fe(t) — (r/i(f))n] = on Fi by (3.3), while 

Wk{t) — Uk{t) = Wh(t) — Uh{t) = on To, integrating by parts we obtain 

(crfe(t) - (Th{t), Ewkit) - Eukit)) = {(Tk{t) - (Th{t),Ewh{t) - Euh{t)) = . 
Therefore, taking the scalar product of both sides of (4.49) with 

{-ak{t),Zk{f))-{-<7h{t),Zh{t)) 
and using (4.48), (4.50), and (4.52) we obtain 

(tTfc(t) - (Th{t), Skit) - eh{t)) + \\Zk{t) - Zh{t)\\l < 

l-t ^ ^ (4.53) 

< ^hk{t){L / Vhk{s) ds + BT{Th + Tk) + C^irr + ^k)} > 
Jo 

where 

i^hkit) := iWcTkit) - cThimi + \\zk{t) - zhimiv^^ ■ 

By (2.29) and (2.30) we obtain tphk{t) < P iphk{t) and 

a^hk{tf < {a-kit) - crh{t),ek{t) - eh{t)) + \\zk{t) - Zh{t)\\l , 
where a := niin{l,ac} and /3 := maxjl, 2/3c} . Therefore (4.53) gives 

a<phk{t) <I3L <phk{s) ds + pB^{Th + Tk) +13C'^{tI'^ + tI'^) . 
Jo 

Using Gronwall's inequality we deduce that 

a^hkit) < l3{B^iTh + Tk) + C^ijl'^ + tI'^)} cMPLt/a) 

for every t G [0, T] and every h and k. As — > by (4.24), recalling (4.51) we obtain that 
Cfe and Zk satisfy the Cauchy condition in L°°{[0,T];L'^{Q;M'^^^)) and L°°{[0,T]; L'^{n)) , 
respectively. □ 

Proof of Theorem 4-4- Step 3. Convergence of the interpolations. By Lemma 4.8 there exist 
two functions e G L^^([0, +00); 2.2(51; Mf^^)) and z G L^^([0, +00); ^^(ll)) such that for 
every T > 

sup ||efe(t)-e(i)||2^0, (4.54) 

0<t<T 

sup \\zk{t) - z{t)\\2 ^ Q . (4.55) 

0<t<T 

Let a{t) := Ce{t) and C(0 := -V'{z{t)). By (2.32) wc have also 

sup \\ak{t) - (T{t)h ^ , (4.56) 

0<t<T 

sup ||C,(i)-C(i)||2-0, (4.57) 

0<t<T 

for every T > . 

By (4.28) and (4.42) for every T > the sequences and arc boimded in 
i?H[0,r];L2(ri;Mf^^)) and H^{[0,T]; L^{n)) , respectively. Hence, up to subsequences, 
we may assume 

e^-e weakly in H^{[0,T];L\n;Mf^^)), (4.58) 
z^ ^z weakly in H\[0, T];L'^{n)) . (4.59) 
Moreover, using (4.42) and the identity 

et{t) = ek{t)+ [ et{s)ds, 



QUASISTATIC EVOLUTION IN PLASTICITY WITH SOFTENING 



23 



which holds for every t € [0, T] , wc obtain 

l|e-(^)-e,(^)||,<4/Vl/Vy^ 

Together with (4.54) and (4.58) this impUcs e(t) = e{t) for every t G [0,r] and that 

sup \\e^{t)-e{t)\\2^0. (4.60) 

0<t<T 

In the same way we deduce that z{t) = z{t) for every t G [0,T] and that 

sup \\z^{t)-z{t)\\2^0. (4.61) 

0<t<T 

This imphes that, without extracting any subsequence, 

e^-e weakly in H^{[0,T];L\fl;Mi^i)), (4.62) 

^ z weakly in H^{[0, T];L'^{n)) , (4.63) 

and that e e H^{[Q,T]; L^{n;M'^,^^)) and z G H\[0,T]; L^{n)) . 

Since div crk(t) = in il and [(Tk(t)n] = on Fi by (3.3), the strong convergence of 
CTkit) to a{t) in L'^{n;M'^^^) yields (ev2)£, taking (2.43) into account. 

By (2.25) and (3.5) we have 

{ptit)'Kit)) = Nlc{<Tk{t)n + {4)n - K~')n,Ck{t) + (1' Ct') a.e. in O 
for <t<t\. Using (4.42) and the Lipschitz continuity of Nf^ we obtain 

{pt it), zt it)) = Nf,{(Tk{t)D, Ckit)) + Rk{t) a.e. in fl , (4.64) 

with ||A/j(f)||2 < C'^T^^'^ , where is a constant, depending on e, but independent of t 
and k. By (4.56) and (4.57) this implies that 

sup \\{pt{t),zt{i))-NI,{<T{t)Dm)h^Q (4.65) 

0<t<T 

for every T > 0. This shows, in particular, that z^ converges strongly in L°°{[0, T]; L^(ri)) 
to a limit that, by (4.61), must coincide with z. Therefore 

sup \\z^(t)-z{t)h^Q (4.66) 

0<t<T 

for every T > 0. Moreover, since (0) = po, (4.65) implies also that there exists p G 
il/„^([0,+oo);i2(O;M^^'')) such that 

sup \\p^{t)-p{t)\\2^0, sup ||p-(f)-p(f) 112^0 (4.67) 

0<t<T 0<t<T 

for every T > , and 

(p(i), z{t)) = N'K{<T{t)DX{t)) a.e. in O 

for a.e. t G [0, +00) . This proves (ev3)e . 

As we did for and e/j , using (4.42) we can prove that 

KW-P.WIl2<cVVvvV2. 

Therefore (4.67) implies 

sup \\p^{t)-p{t)h^O (4.68) 

0<t<T 

for every T > . 

For every t € [0, T] and every k wc have 

Euk{t) = ek{t) + Pk{t) a.e. in 1^, 
Uk{t) = Wk{t) H'^~^-a,.e. in Tq . 
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As Wk{t) w{t) strongly in H^{il;R'^), using the Korn-Poincare inequality, we deduce 
from (4.42), (4.54), and (4.68) that there exists u G i?/„^([0, +oo); //"^(n; M'*)) satisfying 
(evl)e such that 

sup \\ukit) - u{t)\\2 ^ , sup \\Eukit) - Eu{t)\\2 ^ , (4.69) 

0<t<T 0<t<T 

sup \\u^ {t) - u{t)\\2 ^ , sup \\Eu^ (t) - Eu{t)\\2 ^ , (4.70) 

0<t<T 0<t<T 

u^^u weakly in H^{[0,T];H^{n;R'^)), (4.71) 

for every T > . 

Step 4- Uniqueness of the solution. Suppose that {ul,el,pl,zl) and {u"^ , , , z^) are 
two solutions of the £ -regularized evolution problem with boundary datum w and initial 
condition (mo, eo,Po, ^o) • Let cri(t) := Cel{t) , cT'l{t) := 'Cel{t) , Cl{t) := -V'{zl{t)), and 
Ql{t) := -V'izjit)). By (cvS)^ for a.c. t e [0, +oo) wc have 

{Pl{t),zl{t))=NI,{cTl{t)nXl{t)) a.e. inQ, (4.72) 
iplit),zl{t)) = NI,{al{t)D,Cit)) a.e. infi. (4.73) 

Using (evl)£ we get 

{-el{t),zl{t)) = NI,ial{t)D,<:l{t)) - {Eul{t),Q) a.e. inQ, 
{-el{t),zl{t)) = NI,{aUt)D,Clm - {Eul{t),0) a.e. in f2. 

Subtracting term by term we obtain 

i-elit) + el{t),zlit)-zlit)) = 
= Nf,{<Tl{t)n,eM) - NUal{t)nXl{t)) - {Eul{t) - Eul{t),0) ^ > 

a.e. in Q. Since Nf^ is 1/e-Lipschitz and V is M-Lipschitz, there exists a contant L^, 

independent of t, such that 

\\N^^ialit)nXlit))-NI,iaUt)n,Cm\2 < LMelit)-em\2 + Hit)-zUt)h} ■ (4.75) 

Integrating by parts and using the equilibrium condition (ev2)£ and the boundary condition 
iilit) = ul{t) = w{t) H'^~^-a.e. on Fg, we obtain 

{Eul{t)-Eul{t),al{t)-cTl{t))^0. 

Therefore, taking the scalar product of both sides of (4.74) with 

{-al{t) + al{t),zl{t)-zlit)) 

we obtain, using (4.75), 

(alit) - al{t),el{t) - el{t)) + {zl{t) - zl{t),zl{t) - zl{t)) < 
<L,{\\el{t)-em\2 + \\zl{t)-zm\2}\ 

which gives, thanks to (2.29), 

i{Q{elit) - em + hWzlit) - z^U} < UQielit) - + h\\zl{t) - z^ll} , 

for a suitable constant independent of t. 

Since Q(e^(0) - 6^(0)) = and \\zl{0) - z'^^{0)\\?, = 0, from Gronwall's inequality wc 
obtain Q{el{t) - el{f)) = and \\zl{t) - z'^^{t)\\l = for every t G [0, +oo) , hence el{t) = 
eUt), zUt)=zUt), <Ti(t) = o-2(t), and Clit) = Cl{t) for every t e [0,+^). From (4.72) 
and (4.73) we deduce that pl{t) = pl{t) for a.e. t € [0, +00). As p^(0) = pl{Q), this 
implies p\{f) = p1{t) for every t E [0, +00). The equality ul{t) = u^it) follows now from 
the kinematic admissibility condition (evl)£ . □ 
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Remark 4.9. As the solution of the e -regularized evolution problem is unique, the proof 
of Theorem 4.4 shows that it can be approximated by the sequences (wfc, efe,pj., Zfe) and 

(m^ ,e'^,p^,z^) obtained by interpolating the solutions to the incremental problems. The 
precise convergence properties are given by (4.54)-(4.57), (4.60) (4.63), and (4.66) (4.71) . 

Remark 4.10. When d = 1 one can express the conditions of Definition 4.1 in a simpler 
way. 

We begin by observing that condition (cv2)^ is equivalent to the fact that for every 
t G (0, +00) the functions €^{1) and crg(i) are constant on [—5, 5]. 

Let us write A^I-iRxM^MxR componentwise, setting Nf^ = (7V^\7V^^). We select 
an arbitrary representative in the equivalence class of zq. For every y € ["^j^] we can 
consider the solution tp(-,y) of the ordinary differential equation 

dt^{t,y) = N'/{aS),-V'm,y))), ^(0,y) = zo{y) . 

As crs{t) and V'{'ip{t, y)) are bounded uniformly with respect to t and y, the same property 
holds for dtil>{t, y) . By the dominated convergence theorem we have 

^{t + h.-) - r{l.-) 

> dti:{t,-) 

strongly in L'^{\—\,\]), so that t ■0(f, •) is a solution of the Cauchy problem in 

m = N'/{<Te{t), V'(V(i))) . V'(O) = zo . 
By condition (evS)^ the function is a solution of the same Cauchy problem. Since 
the function xjj N^'^{(Tir{t),iJj) is a Lipschitz continuous function from L^([— i,^]) into 
I/^([— ^, ^]), by uniqueness we conclude that for every f > we have Ze{t,y) = tp{t,y) for 
a.e. y G [-i, i] . 

This shows that for every t > we can choose representatives of Zg{t) and z^{t) in their 
equivalence classes such that z^{t,y) — dtZs{t,y) for every y G [— i, i] . Using (evl)e, for 
these representatives condition (ev3)e is equivalent to 

ee{t) = - j\ ^K'K(t), Ce(i,2/)) dy + w{t, \) - w{t,-\) , (4.76) 

zS, y) = N'/{(tS), CAt, y)) for every y g (-§, |) . (4.77) 

If, in addition, zq G C'([— 55!]), then z^ G C([0, +00); C([— i, i])) , since (Te is continuous 
and z^{-,y) is a solution of an ordinary differential equation with initial condition Zo{y) 
depending continuously on y . 

5. Approximable quasistatic evolution 

In this section we state and prove the main results of the paper on approximable qua- 
sistatic evolutions. 

5.1. Definition and properties. We recall that a function v from [0, +00) into the dual 
F of a Banach space is weakly* left-continuous if 

v{s) v{t) weakly* in Y 

as s ^t, with s <t. We recall that a system /j, G 5Gy([0, -|-(X)),rT;M^^'^xR) is said to be 
weakly* left- continuous if for every finite sequence ti, . . . ,tm in [0, -|-oo) with ti < ■ ■ ■ < tm 
the following continuity property holds: 

At.,....„ ^ f^t,...t^ weakly* in GY{n; (M^><'^xM)") (5.1) 

as Si ^ ti, with Si G [0, +00) and Si <ti. 
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The notion of approximable quasistatic evolution is made precise by the fohowing defini- 
tion. 

Definition 5.1. Let w e Hl^{[Q,+oo);H^{9.;M.'^)) , uq e H^{0), eo e L'^{Q.;Mi^^) , 

Po G M^^'*) , and zq G L'^{Q). Suppose that the kinematic admissibihty condition 

(evl)o and the equihbrium condition (ev2)o of Theorem 4.4 are satisfied. An approximable 
quasistatic evolution with boundary datum w and initial condition {uq, eo,po,zo) is a triple 
(w,e,/x), with 

u e Lf^^{[0,+^y,BD{n)) , e e L,- ([0,+oo);i2(Q;M^;^)) , 
SGY{[0,+oo),Ti; M^^'^xR), 

and u, fi weakly* left-continuous, for which there exist a positive sequence cfe ^ and 
a set 6 C [0, -|-oo), containing and with £^([0, -|-oo) \ 0) = 0, such that the solutions 
(ttg^ , e^^ , Pgj, , ^ej. ) of the £fe -regularized evolution problems with boundary datum w and 
initial condition {uo,eo,po, Zq) satisfy 

{t) u{t) weakly* in BD{n) , (5.2) 

e,, it) - e{t) weakly in L^{n; Mf^^) (5.3) 

for every t G @, and 

(<^(p.,,-e,))t,...*^ ^ Mt,...t. weakly* in GY{n; (M^x-^xR)™) (5.4) 

for every finite sequence ti, . . . ,tm in © with ti < ■ ■ ■ < tm ■ 

If, in addition, fi = <5(p_z) for some functions p G L^^([0, -hoo); M,,(r2; M^'''')) and z G 
-^toc(P' -^^(^)) ; we say that (tt, e,p, z) is a spatially regular approximable quasistatic 

evolution. 

The first estimates for approximable quasistatic evolutions are proved in the next lemma. 

Lemma 5.2. Let w , uq, eo, Po, zq, (u,e,/Lt), £k, and {ue^,ee^,p^i^, z^^) be as in Defi- 
nition 5.1. Then for every T > there exists a constant Ct < +oo such that 

sup \\ee,{t)h<CT, VH{p,^,Ze,;0,T)<CT, (5.5) 

«e[o,T] 

£fc / \\PeM\2dt<CT, and Ek [ \\z,,{t)\\ldt<CT. (5.6) 
Jo Jo 

Proof. By the energy equality (ev4)£ and by (4.22) for every T > we have 
Q(e,, (T)) + Vnipe, , z,, ; 0, T) + V(z,, (T)) 

[ \\p,M\ldt + £kl \\zeM\\ldt= (5.7) 
Jo Jo 

= Q{eo)+V{zo)+ [ {cT,,{t),Ew{t))dt, 
Jo 

where cr^,(t) —Ce^Ji). 

The proof can now be concluded as in Step 2 of the proof of Theorem 4.4. □ 

To study the properties of approximable quasistatic evolutions we need to introduce some 
definitions. Given /j, G SGY{[0,+oo),Tl;M'j^'^xR) , its dissipation 'DH{lJi;a,h) on the time 
interval [a, b\ C [0, -|-oo) is defined as 

k 

sup^(iJ(^i - ^i-i,6i - 9i-i), Htou...tki^, Co, do,..., Cfe, ek,v)) , (5-8) 

i=l 
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where the supremum is taken over ah finite famihes to,ti, . . . ,tk such that a = to < ti < 
■ ■ ■ <tk = h. As in the case of the variation, we have 

k 

VH{n\a,h) = sup^(ff(^i -^i_i,6li - ei-l),|x^._^^^{x,^,i-l,ei-l,£,i,0i,r|)) , (5.9) 

i=l 

where the supremum is taken over all finite families to,ti, . . . ,tk such that a = to < ti < 
■■■ <tk = b. _ _ 

li fj, = 5(p,2) for some p: [0,+oo) ^ Mb{n;M'j^'') and z: [0, +00) ^ Mb{n) , we have 

VH{ti;a,b)=VH{p,z;a,b), (5.10) 

where 'Dh{p, z;a,b) is defined in (4.21). 

We also introduce the homogeneous functions {V}: KxM — > M and {F'}: MxR — > R 
defined by 

\v{e/ri) if?7>0, 
F°°(6') if??<0, 



and 

{V'}{6,r,): 



vV'{0/r]) ifr]>0, 

if ?7 < . 

According to [6, Definition 2.15], for every fi e GY{n;M'j^'^xR) we can define the 
measure T^Q{{V'}fj,) e Mb{Q) as 

(<^, 7rj^({y'}/x)) = (^(a;){y'}(^, 77), /x(a;, ^, 6, rj)) (5.11) 

for every (p G C(f2) . 

Lemma 5.3. For every ^ e GY{Q;Mj^'^xR) we have ir^UV'jfj,) e i°°(0) and, moreover, 

lkn({nM)lloo<!|r||eo. 

Proof. By (2.1) we have 

|((p(a;){n(^,^),M(-^,e,^^,^))l < ||V^'||co(|^(x)|r,,M(x,e,e,r/)) = ||l^'||oo||^||i 

for every ip € C(0). The conclusion follows. □ 

The main properties of approximable quasistatic evolutions are proved in the following 
theorem, where )C{fl) is the convex set defined in (4.5). 

Theorem 5.4. Let w, uo, eo, Po, o,nd zq be as in Definition 5.1. Let {u,e,fj,) be an 

approximable quasistatic evolution with boundary datum w and initial condition {uo, eo, /Uo) > 
and let 

(p(t),z(t))=bar(/x,), a{t):=Ce{t), C{t) := -n^{{V'}Ht) . (5.12) 
Then the following conditions are satisfied: 

(evl) weak kinematic admissibility: for every t £ (0, +00) 

Eu{t) = e{t) + p{t) in U , 
p{t) = {w{t) - u{t)) © nH'^-i on To ; 

(ev2) equilibrium condition: for every t G (0, +00) 

div<T(f) =0 in fl, [f(t)n] =0 on Fi ; 



(ev3) stress constraint: for every t E (0, +00) 

{cTn{t),<:{t))&m); 
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(ev4) energy inequality: for every T € (0, +00) 

Q(e(T)) + 0, T) + {{V}{e, 7?), 9, r,)) < 

< Q(eo) + V{zo) + I {(T{t), Ew{t)) dt . 

Jo 

Remark 5.5. If (u, e,p, z) is a spatially regular approximable quasistatic evolution, then 
(^{t) = —V'{z{t)) by (2.6), (5.11), and (5.12). Therefore the energy inequality (ev4) reduces 
to 

Q(e(T)) + V„{p, z; 0, T) + V(z(T)) < 

(5.13) 

< Q(eo)+V(^o)+ / {<T{t),Ew{t))dt 

Jo 

thanks to (5.10). 

Proof of Theorem 5.4- Let Q, Sk, and (ttg^ , , , ) be as in Definition 5.1. By (2.5) 
and (5.2)-(5.4) we have 

Wefe(i) ^ u{t) weakly* in BD{Q.) , 
ee,{t) - e{t) weakly in (fi; Mf^>=^) , 
p^^{t) p{t) weakly* in Mb{Ti; M^^**) 

for every t G 0. By applying [5, Lemma 2.1] wo prove that (evl)^ imples (evl) for every 
t £ Q. These equalities are then extended to every t G (0, +00) by left-continuity. 

For every t G Q the equilibrium condition (ev2) follows from the equilibrium condition 
(cv2)£ satisfied by (Tg^{t) and from (2.43). The same result can be obtained for every 
t G (0, +00) by left-continuity. 

By the modified stress constraint (evS)^ , for every k and a.e. t G (0, +00) we have that 

icT,,it)D - ekpJt)Xe,{t) - SkZeM G IC{n) . (5.14) 

For every t € Q we have that J(p^ (f),ze^.(t)) ^ A*t weakly* in Gy(0; M^^'^xM) , which, by 
(2.5) and Lemma 5.3, implies that Ceji^) ~^ weakly* in L°°{Q). By (5.6) the sequences 
^fc||P£fc(')l|2 and efe||i;efc(-)||2 converge to strongly in _L^^^([0, -|-oo)) , hence there exists a 
sequence kj oo such that SkjPe^.{t) ^ in L'^{n;M'^'^) and ekjZe^,{t) ^ in L'^{Q) 
for a.e. t G [0, -|-oo) . This implies 

{a{t),C{t))elC{n) (5.15) 

for a.e. t G (0, -|-oo) . 

The left continuity of /i, together with Lemma 5.3, implies that the function (J is weakly* 
left-continuous in L°°{n). As cr is left-continuous too, by approximation (5.15) holds for 
every t G (0, -|-oo). This concludes the proof of (ev3). 

It remains to prove the energy inequality (ev4). By (5.3) and by the lower semicontinuity 
of Q , for every T G 6 we have 

Q(e(T))< hminf Q(ee,(r)), (5.16) 

k^oo 

and, by left-continuity, for every T G [0, +00) we have 

Q(e(T))<liminfQ(e(5)). (5.17) 

S<T 

As the dissipation is lower semicontinuous (see [6, Theorem 8.11]), using (5.4) and the 
left-continuity, and taking into account (5.10), one can check that for every T G 

VH{f^;0,T)< hminf %(p,^,z,,;0,T), (5.18) 

fe— >(X) 
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and for every T G (0, +00) 



Vh{ij.;0,T) <\immiVH{fJ.;0,S) . (5.19) 

S<T 



Moreover, for every T e 6 we have 

V(ze, (T)) ^ {{Vm v), I^Ax, ^, V, 6)) (5.20) 

as fc — > 00, and for every T £ (0, +00) 

{{V}{e, 77), txsix, V, 6)) - {{V}{9, v), t^A^, ^, V, 0)) (5.21) 

as 5->T, with S<T. 
Finahy, by (5.3) we obtain 

{rT,Xt),Ew{t)) ^ {(T{t),Ew{t)) 

for a.e. t G [0, +00) . Thanks to (5.5) we can apply the dominated convergence theorem and 
we obtain 

/ {creJ(b),Ew{t))dt^ [ {cr{t),Ew{t))dt. (5.22) 

Combining (5.16), (5.18), (5.20), and (5.22), we can pass to the limit in the energy equality 

(ev4)e satisfied by (m^^ , e^^ , p^^ , z^^ ) and we obtain the energy inequaUty (cv4) for every 
T G 6 . The result for every T G [0, +00) can be obtained by left-continuity thanks to 
(5.17), (5.19), and (5.21). □ 

5.2. The existence result. We now prove the main existence result of the paper. 

Theorem 5.6. Let w G iJ/„^([0, +00); iJi(^^; K'^)) . G H^{n), eo G L'^{9.;Mi^^) , 
Po G L^(ri; M^^'^) , and zq G L^(ri). Suppose that the kinematic admissibility condition 
(evl)o and the equilibrium condition (ev2)o of Theorem 4-4 f^fc- satisfied. Then, there ex- 
ists an approximable quasistatic evolution with boundary datum w and initial condition 
{uo,eQ,pQ,Zo) . 

Proof. Let us fix a positive sequence — > 0. For every k let (tie^ , ee,, , p^^, , Ze/t ) be the 
solution of the 6^ -regularized evolution problem with boundary datum w and initial con- 
dition (uq, eo,po, Zq) , and let crg^{t) := Cse^it). From the energy equality (ev4)£ for every 
T > we obtain 

Qie,,iT)) + VHip,,,z,,;0,T) + V{z,,) < Q{eo) + V{zo) + I {a,,{t),Ew{t)) dt . 

Jo 

By (2.11) and (2.35) we have 

Vh{p,„z,,;0,T) + V{z,,) - Vizo) > 
> C^Var(p,^,z,,;0,T) = C^Var(5(p^^,,^^); 0, T) . ^^"^-^^ 

Therefore, 



Q(ee.(T)) < Q(eo)+ Ms), Ew{s)) ds . 
Jo 



Arguing as in the proof of (4.41) and (4.42), we obtain that there exists a constant Ct, 
depending on T but independent of t and k , such that 

||es.(t)||2 <Ct, M)h<CT, (5.24) 
Var (<5(p^^ ) ; 0, T) = Var (p,^ , z,, ; 0, T) < Ct , (5.25) 

for every t G [0, T] and for every k. Inequality (5.25), together with the initial conditions, 
implies that ||p£^(t)||i < |bo || 1 + ||-2o||i + C't and ||z£j^(i)||i < ||po || 1 + ||-2o||i + C't for every 

tG [o,r]. 
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Using (2.12) and the kinematic admissibility condition (evl)e of Definition 4.1, from the 
estimate of ||Pe^(i)||i and from (5.24) we deduce that for every T > the sequence Ue^{t) 
is bounded in BD{^1) uniformly with respect to t E [0,T] and k. 

By applying Helly's theorem for systems of generalized Young measures (see [5, Theo- 
rem 8.10], together with a standard diagonal argument, we construct a subsequence, still 
denoted Sk, & set 9 C [0, +oo) , containing and with [0, T] \ 8 at most countable, and a 
left continuous n € 5'Gy([0, +oo),l7; M^'^'^xR) , with 

Var(/i; 0, T) < Ct , (5.26) 

llMtll* < lbo||i + |No||i+CT for every t e [0,T], (5.27) 

such that 

(<5(p,^,..j)t....t^ ^ Mt,...t„ weakly* in GF(0; (M^x'^xM)'") (5.28) 

for every finite sequence ti,. . . ,tm in © with ti < ■ ■ ■ < tm- For every t G [0, +oo) we set 
{p{t),z{t)) := bar(/xj. 

Let us fix i e 6 . Since u^^ (t) and e^^ (t) are bounded in BD{U) and ^^(O; Mf^^) , there 
exist an increasing sequence kj (possibly depending on t) and two functions u{t) € BD{fl) 
and e{t) e L2(ri;Mf^„^J such that u^^^{t) u(t) weakly* in SD(fi) and e^^^^) e{t) 

weakly in L^(f2; M^^^) . As in the proof of Theorem 5.4 we can show that u{t) , e{t) , p{t) , 
and w{t) satisfy the weak kinematic admissibility condition (evl). 
By the equilibrium condition (ev2)e we have 

div (Tsjl^t) = in O and [f£fc(^)?^] = on Fi . 

Using the weak definition of divergence and normal trace (2.43), we can pass to the limit 
and we obtain 

div a{t) = in O and [<T{t)n] = on Ti . (5.29) 
For every ip e H'^{n;M^) with ip = H^-'^-a.e. on Tq we have 

Q{e{t) + E^) - Q{e{t)) = {a(t), E^) + Q{E^) = Q{E^) , 
where the last equality follows from (2.43) and (5.29). Therefore 

Q(e(i)) < Q{e{t) + Eip) . (5.30) 

If kj is another sequence such that u^. (t) u(t) weakly* in BD(Cl) and e^. (t) 

e{t) weakly in L^{Cl;Mfy^) , then u{t), e{t), p{t) , and w{t) satisfy the weak kinematic 
admissibility condition (evl) and 

Q{e{t)) < Q{e{t) - Eip) (5.31) 

for every ip G H'^{Q:R'') with ip = W^-^-a.e. on Tq. 

From the weak kinematic admissibility condition (evl) it follows that u{t) — u{t) = 
W'-i-a.e. on Tq, and e{t)-e{t) = Eu{t)-Eu{t) a.e. on O, hence u(t)-u{t) e H'^{Q.;W^). 
If we take (p = ^{u{t) — u{t)) in (5.30) and (5.31), by adding the inequalities we obtain 

Q(e(i)) + Q(e(0)<2Q(i(e(t)+e(t))), 

which implies that e{t) = e{t) a.e. on by the strict convexity of Q. Therefore u{t) = u{t) 
a.e. on fi, since Eu{t) = Eu{t) a.e. on Q and u{t) = u{t) H'^~^-&.e. on Fq. As these limits 
do not depend on the subsequences, the convergence results hold for the whole sequences. 
This proves (5.2) and (5.3). Moreover, (5.24) and the uniform estimate on UeJ(t) in BD{Cl) 
give 

sup ||e(t)||2 < +0O and sup [[^^^(i)!!! < +oo (5.32) 

teen[o,T] teen[o,T] 

for every T > 
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We now show that for every t G we have 

(s) ^ e(t) weakl; 
u{s) — ^ u{t) weakly* in BD{fl) 



e(s)-e(t) weakly in L2(f^;Mf^^), ^^^^^^ 



symJ 



as s ^ t, with s S O and s <t. By (5.32), for every t ^ Q there exist e e L (f2 
u e BD{Q) , and a sequence Sj G 0, with sj t and Sj < t, such that e(sj) — ^ e weakly 
in L^(f2;Mg^^) and u{sj) u weakly* in BD{fl). As each e{sj) satisfies the minimality 
property (5.30), we deduce 

Q(e) < Q(e - Eip) (5.34) 
for every (p G ff-^(r2;R'') with = H'''~^-a..e. on Fq. Moreover by left-continuity we have 
p{sj) p{t) weakly* in Mi>,(Q; M^^*^) , which implies that w, e, p(f), and w{t) satisfy 
the weak kinematic admissibility condition (cvl) of Theorem 5.4. As before wc can take 
(fi := ^{u — u{t)) as a test function in (5.30) and (5.31) to deduce from the strict convexity 
of Q that e = e{t) and, in turn, « = u{t). Since the limit is independent of the sequence 
Sj , wc have proved (5.33). 

The same argument shows that for every t € [0,T] we can define in a unique way e(t) € 
L2(f2;Mf^^) and u{t) e BD{n) such that (5.33) holds as s ^ i, with s e [0,T] and 
s<t. □ 

Remark 5.7. Assume that w, Tq, uq, cq, po, zq satisfy (4.12) (4.16). Then by Proposi- 
tion 4.5 the solutions of the £-regularized evolution problems are spatially regular. There- 
fore every approximable quasistatic evolution with boundary datum w and initial condition 
{uo, eo,po,zo) is spatially regular and has the special form 

uit,x) = mx, eit,x) = e{t), p{t,x) = e{t), z{t,x) = e{t), (5.35) 

with ^« : [0, -1-00) ^ Mf^4 , : [0, +00) M^'''' , 61 : [0, -|-oo) ^ R left-continuous and 

c{t) + e{t)^e{t) (5.36) 

for every t € [0, -|-oo) , where ^^{t) is the symmetric part of ^{t) . 

6. The case of simple shear 

We analyze in this section the case of simple shear in the Dirichlet-Periodic case (DP) for 
an isotropic material in dimension d > 2 with shear modulus /x (see (2.27)). We consider 
^-periodic solutions with boundary data of the form 

w{t,xi,x) := \/2w^{t,xi)e2 (6.1) 

and initial conditions of the form 

uo{0,xi,x) := \/2?i^(xi)e2 , eo{0,xi,x) := M(e^) , 
Po{Q,xi,x) := M{p^{xi)) , zo{Q,xi,x) := z^{xi) , 

where w'^ € i?/,,([0, +00); i?i([-i, i])), € H\[-\, i]), eR,p^, z^ G L^{[-^, i]), 
62 is the second element of the canonical basis of M.^ , and M:'M.—^ M.j^'^ is the linear 
isometry defined by 

/ ^ ••• o\ 

••• 



(6.2) 



M{a) := 



V2 



(6.3) 



\ ••• / 

which connects the amount of a pure shear with its strain. 

We assume that the elastic domain K depends on the stress a only through its euclidean 
norm \a\. More precisely, we assume that there exisits a closed convex set c MxM such 
that 

K = {(a,C) G M^-'^xR : (|a|,C) G K^} ; (6.4) 



32 



G. DAL MASO, A. DESIMONE, M.G. MORA, AND M. MORINI 



it is not restrictive to assume that 

{a,OeK''' ^ (_c,C)ei^^'. (6.5) 

Theorem 6.1. Assume (6.1), (6.2), (6.4), and (6.5). Then (Ug , Cg , , ) is the solu- 
tion of the e -regularized evolution problem with boundary datum w and initial condition 
{uo, eQ,pQ,zo) if and only if 

Ue{t,xi,x) = V2u^{t,xi)e2, e,{t,xi,x) = Mie^it)), (6.6) 
p,{t,xi,x) = M{p^{t,xi)), Ze{t,xi,x)^z^it,xi), (6.7) 

where [u^J .ef ,pf , z^) is the solution of the e -regularized evolution problem, in dimen- 
sion d = 1, corresponding to = with boundary datum and initial condition 
{uq", Cq-jPo", Zq) . 

Proof. Let us first prove that for every (a, C) € KxR we have 

PK(M(a),C) = (M(d),C) ^ PKn{a,0 = {a,0. (6.8) 

Indeed on the one hand, if (q;,C) := -Pft:«(ajC): then the distance between (M(a),^) and 
{M{a),C,) is larger than the distance between {M{a).Q and {M{a)X). This is an easy 
consequence of the fact that Af is a hnear isometry. On the other hand, if (a, C,) := 
PK{M{a),() , it follows from (6.4) that a is the projection of M{a) onto the ball {a £ 
M'^'' : \a\ < \a\}, which implies that a = M{a) for some a e M. These two facts together 
allow to establish (6.8). 

Using the linearity of M we can prove that (6.8) implies 

(M(d),C) = iV^(M(a),C) ^ (d,C)=iV^H(a,C) (6.9) 

for every a , , d , C in K • 

By (6.2) the kinematic admissibility condition (evl)o (Theorem 4.4) for {uo,eo,Po,Zo) 
is equivalent to (cvl)o in dimension d = 1 for (u^, e^, p^, Z(f ) , while the equilibrium 
condition (ev2)o is always satisfied, since cq and Cq are constant. Assuming (evl)o, let 
{uf,ef,pf,zf) be the solution of the one dimensional e-regularized evolution problem 
with boundary datum and initial condition [u^ ,p^ ^ z^) ^ and let (ttg, e^, p^, ^g) be 
defined by (6.6) and (6.7). It is easy to check that {ue,ee,p^, Zg) satisfies the initial condi- 
tions (evO)^ of Definition 4.1. Moreover, the kinematic admissibility (evl)^ follows from the 
fact that E{u^(t)) = M{Duf (t)) , and from the kinematic admissibility in dimension one 
{D denotes the distributional derivative with respect to the space variable). Equilibrium 
condition (ev2)£ follows from the fact that (T^{t) is constant on (— ^, ^) by Remark 4.10. 
Finally the regularized flow rule (evS)^ follows from the condition in dimension one and 
from (6.9). □ 

Remark 6.2. The onc-dimcnsional problem has also different interpretations, for example 
in the study of uniaxial loading of cylindrical bodies. In this case, the symmetry condition 
(6.5) is no longer a natural assumption since one may envisage a different behaviour in 
tension and compression. Theorem 4.4 holds true also in this general case, since the only 

hypotheses on K arc (2.15) and (2.16). 

In order to prove the equivalence between the approximable quasistatic evolutions for 
simple shears and the solutions to the corresponding reduced one-dimesional problems, for 
every m we consider the linear operator iprn- (RxM)™ ^ (M^^'^xR)™ defined by 

V'm((/?l, ^l), • • • , {I3m, 0m)) = {{MiPl), 61), {M{pm), 0m)) ■ 

In the following theorem wc shall prove that the system of Young measures fi in the ap- 
proximable quasistatic evolution is related to the system Young measure /j,^ of the reduced 
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problem by the formula 

{f,IJ't,...tJ =j{f{{xi,x),iJrn{i),r]),ixf^_t^{xi,i,ri))dx (6.10) 

for every / G C''°"(Qx (MxM^'''')™ xM) and every finite sequence ti,...,tm in [0, +00) 

with ti < ■ ■■ <tm- 

Theorem 6.3. Assume (6.1), (6.2), (6.4), and (6.5). Then {u,e,fi) is an approximable 

quasistatic evolution with boundary datum w and initial conditions (uq, eo,Po, Zq) if and 
only if (6.6) and (6.10) hold and {u^,e^,iJ,^) is an approximable quasistatic evolution for 
the problem in dimension d = 1. corresponding to = with boundary datum and 
initial condition (uq , Cq , p^ , Zq ) . 

Proof. By Theorem 6.1, the solution {us,ee,Pg, z^) of the e-regularized evolution problem 
with boundary datum w and initial condition {uq , cq , Po , zq) satisfies (6.7). Therefore, if 
{u^, e^,p^, zf ) is the solution of the reduced e-regularizcd evolution problem with bound- 
ary datum and initial condition (uqjBqjPqjZq), then (6.10) holds with fj, replaced by 
^{p^,ze) 9'nd replaced by ^(pH ^H). The conclusion follows. □ 

7. A SPATIALLY HOMOGENEOUS EXAMPLE 

In this section we assume that d > 2 and C is isotropic, which implies that 

for some constants fj. > and k > 0. We also assume that 

K := {(^, e) e M^^^-^xM : + 6»2 < 1} , (7.1) 

V{0):=l-lVT+¥, To-dfl, ri:=0. 

Let us fix a constant 9o > and a dxd matrix ^0 with tr^o = 0. Wc assume that the 
symmetric part of ^0 is different from . We will examine the approximable quasistatic 
evolution corresponding to the boundary datum 

w{t, x) := t^ox , 

and to the initial conditions 

uo{x)=0, eo(a;)=0, po{x)=0, zo{x) = do- 

7.1. The special form of the solution. The following theorem shows that, in this case, 
any approximable quasistatic evolution is spatially regular and satisfies u{t) = w(t) for 
every t G [0, +00), while e{t) , p{t) , and z{t) do not depend on x. Moreover, e{t) and p{t) 
are proportional to the symmetric part Q of ^0 • 

Theorem 7.1. Let C, K, V, Fq, Fi , Co; Co; ? "O; cq, po, and zq satisfy the 
assumptions considered at the beginning of this section. Then any approximable quasistatic 
evolution with boundary datum w and initial condition {uo,eo,po,zo) is spatially regular 
and has the special form 

u{t,x):=t^ox, e{t,x) ■.= {t-ij{t))eo, p{t, x) := , z{t,x)=e{t), (7.2) 

where tp: [0, +00) [0, +00) and 6: [0, +00) [0, +00) are functions depending on fi and 
60 . We have ■tp{t) = and 0{t) = 9o for <t < to, where 



while 

1 



V'(0=i-^\/l-^'(^(i))' (7.3) 



for t> to. 
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Proof. By Proposition 4.5 the solution {ue,ee,Pg,Ze) of the e -regularized evolution problem 
with boundary datum w and initial condition {uQ,eo,po, zq) has the form 

u,{t, x) = w{t, x) , eS, := {t - Ve(t))Co . ^) := Mt)^o > ^S, x) := Beit) , 

where ipg, e ^/^^([O, +00)) is the unique solution of the Cauchy problem 

{Mmo,m) e Ni,{2f,{t-Mtmo,-v'ie,m (7.4) 

V'e(o) = o, e,{o) = eo. (7.5) 

According to (2.24) and (7.1) we have 

iV^(a,C) = i^%£=(v^^^- 1)+. (7.6) 

As Nj^ is Lipschitz continuous we have ips, Os & C^([0, +00)) . It is easy to see that in the 
interval [0,to\ the pair (^^(t), ^^(t)) := (0, ^0) is a solution to (7.4) and (7.5). 
By Remark 4.2 we have 

{2^i{t - - eMmo, -V'{9,{t)) - s9,{t)) G dH{Mtn, Os{t)) (7.7) 

for every t G [0, +00) . Let us prove that 

t-^e{t)>0 (7.8) 

for every t > to. We argue by contradiction. If (7.8) is not true, let 

ti := M{t >tQ:t- ip{t) < 0} . 

Since V'e(^o) = 0, by continuity we have ti > to, ti — tljg{ti) =0, and t — ^e{t) > for 
t e [to,ti). It follows that 1 - ipe{ti) < 0. On the other hand, (7.4) and (7.6) imply that 
"fpeiti) =0, which contradicts the previous inequality. This concludes the proof of (7.8). 
We now show that for every t > to we have 

{2^x{t - - sM^o, -V\0S)) - ee,{t)) € dK , (7.9) 

which is equivalent to 

\C^n2fi{t - Mt)) - sMt)? + [v'ieM) + = 1 . (7.10) 

Arguing by contradiction, let us consider a maximal open interval {ti,T2) in which (7.9) 
does not hold. Since dH{^, 6) C dK for every (^, 6) ^ (0, 0) , we deduce from (7.7) that 

tps{t) = Oeit) = for every t d (ti, r2) . Therefore there exist two constants Ci and C2 such 
that il^eit) = ci and 9^{t) = C2 for every t £ (ti,T2). It follows that 

{2fi{t - i;,{t))eo - ei^,{tn,-V'{e,{t)) " eUt)) = {2t,{t - ci)^^, -F' (ca)) , 
which by (7.7) gives 



/4n-^{t - Ci)2|^o-|2 + T/'(c2)2 < 1 (7.11) 
By the maximality of {ti,T2) we have (2/x(ti — ci)^q, —V'{c2)) & dK and hence 

y'4/i2(ri-Ci)2|e^|2 + y'(c2)2 = l. 

On the other hand, on the interval (ti,T2) we have t — ci > by (7.8). This implies that 

^4M2(i-Cl)2|^^|2 + y'(c2)2>l 

for every t G {ti,T2) which contradicts (7.11) and concludes the proof of (7.9). 

Using Definition 5.1 and the special form of the function {u£,ee,Pg,Ze) , it is easy to 
prove that there exist two left-continuous functions ip: [0, +00) — > M and 6: [0,-|-oo) — > M 
such that 
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for a.e. t G [0, +00). This property, together with the weak* left-continuity of e, p, and 
z, implies that (7.2) holds for every t > 0. Using (5.6) we obtain also ekipsk{t) ~^ and 
efe4.(i)->0 for a.e. t>0. 

Passing to the limit in (7.10) and using again left-continuity, we get 

for every t > 0. Since t - V'(i) > by (7.8), we obtain (7.3). □ 



7.2. A second order equation for the regularized evolution. Thanks to (7.2) and 
(7.3) the behaviour of an approximable quasistatic evolution is completely determined by a 

function 9 , which is the limit of a sequence of funtions 6^ related to the e -regularized evo- 
lution problem. The following lemma provides a careful analysis of a second order equation 
satisfied by the functions . The special properties of the different terms appearing in this 
equation will be important in the proof of the asymptotic behaviour of as e ^ 0, and 
consequently in the description of the limit function 6 . 

Lemma 7.2. Let £ > and let 9^ he the function introduced in the proof of Theorem 1.1. 
Then 9^ is the solution of the Cauchy problem 

e9e = F,{9„9,) on [to, +00), 9^0) 4(io)=0, (7.12) 

where 

Fe{9, 9) := A{9, e9) + B{9, e9)9 , (7.13) 

A{9,v) := Ao(^)- Al(^,^;)^; + A2(^,^;)^;^ (7.14) 

A^{9) := -2p.\a\V\9),Jl-V'{9)\ (7.15) 



A^{9,v) := M^^Wl-{V'i9)+v) 



+ 2.|^g|^-W ^^-^^^^^^^t-^^^^ , (7.16) 



A2i9,v) := -^^l-{V'{9)+v)\ (7.17) 
B{9, v) := -Bo{9) + Bi{9)v + B2{9y - B3{9y , (7.18) 



Bq{9) := 2/z(l - V'{9y) + V'{9yv"{9) , (7.19) 

Bi{e) .- yr^^ , (7.20) 

B2{9) :=3{2fi-V"(9)), (7.21) 
_ 2n-V''i9) 

- — W{9r~ ■ ^ ' 

Moreover, 9s{t) > for every t > to and 0e{t) +00 as t ^ +00. 

Remark 7.3. We observe that for every 9 > 9o and every v G [0, — ^'(6*)] all functions 

Aq, Ai, A2, Bi, B2, B3 are positive, except Bq. In particular, in this range the following 
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inequalities hold: 

-x/3m|Co1^'(^^o) < A„{e) < ^i\eo\ , (7.23) 
2^,ro\<A,{e,v)<4^i\^^o\, (7-24) 

2V3f,\eo\<A2ie,v)<-p^^, (7.25) 

lfi-l<Bo{e)<2fi, (7.26) 

6H < B2{0) < 6ii + ^ , (7.28) 
il^<B3{9)<-yj^{2,, + ^). (7.29) 

Proof of Lemma 7.2. Let V-'e be defined as in the proof of Theorem 7.1. We observe that 
(7.4), (7.6), and (7.8) imply that V-'e(i) > for every t > to ■ Let us prove that also 

6e{t)>0 for every t > io • (7.30) 

If not, we can define t2 = infi;, where E := {t > to : 9e{t) < 0}. By (7.4) and (7.6) we 
have —V'{6s{t)) < for every t G E. By continuity, 

-V'{e,{t2))<0. (7.31) 

On the other hand 6^{t) > for every t G [to,t2]. Hence < 6q = 9^{tQ) < 9s{t2). As 
—V is increasing, we dedude that < —V'{9q) < V'{9s{t2)), which contradicts (7.31) and 
concludes the proof of (7.30). 

By (7.9), the fiow rule (7.7) is equivalent to saying that for every t > to there exists 
X{t) > such that 

A(t)(2/i(t - ^e(t)) - sMt)) = Mt) , (7.32) 
X{t){-V'{9S)) - £0e{t)) = 9S) ■ (7.33) 
To simplify the notation it is convenient to introduce the function : M.xM. M. defined 

by 

W,i9,y) ■.= V'{9)+ey. 
We shall frequently use the inequalities 

-^<W,{9,y)<0 

for every 9 > and y e [0, —V'{9)/e) . Let us prove that 

2fi{t - i;,{t)) - eij,{t) > , WM{t),&s{t)) <0 (7.34) 

for every t>to- H ipe{t) > 0, the first inequahty in (7.34) follows from (7.32). If tps{t) = 0, 
the same inequality follows from (7.8). If 6g{t) > 0, the second inequality in (7.34) follows 
from (7.33). Finally, as 9^ is nondecreasing by (7.30), we have < -V'{9o) < -V'{9s{t)) 
for every t>to- This implies the second inequality in (7.34) when 0e{t) = 0. 
Prom (7.33) and (7.34) we obtain 

xit) = . 

W,{9e{t),9,{t)) 

Substituting X{t) in (7.32), we have 



W,i9S),0S)) ^'V W,i9S),0e{t))^ ' 'We{9S),0S)) 
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By (7.10), (7.34), and (7.35) we obtain 



1 J I -W^ {9, {t), Kit)) . 

l?ol We{9e{t),ee{t)) 



From (7.35) it follows that 



On the set where 9e{t) we obtain 



2,(,-^.W) = '-''!'f»V'-^''-'^'"'-'''"». (7.37) 



On the set where = (7.37) reduces to 



2fi{t - Mt)) = • 

which is a consequence of (7.10) and (7.34), taking into account that ■^e(t) = where 
9s{t) — (this can be easily deduced from (7.32), (7.33), and (7.34)). 

By differentiating (7.37) and using (7.36) we finally obtain that 9^ satisfies the equation 
in (7.12). 

Since 9^ is of class on [0, +00) and 9e{t) = for t e [0, to) , we deduce that 9e{to) = 
and 4(to) = 0. 
Let us prove that 

9s{t) > for every t > to ■ (7.38) 

First we observe that the equation gives 9e{to) > 0, so 9e{t) > in a right neighbourhood 
of to • If this inequality is not satisfied for every t > to, let r be the first time greater than 
to such that ^^(t) = 0. From the equation we get 9^{t) = F^{9e{T),0) = Ao{9e{T)) > 0, 
which leads to a contradiction with the definition of r, and concludes the proof of (7.38). 
Let 

9:= lim 9eit) . 

We want to prove that 9 = +00. If not, by (7.34) and by the monotonicity of 9^ , we have 
that 

< e9e{t) < -V'{9e{t)) < -V'{9) , ^0 < 9,{t) < 9. 
This implies that the pair {9e{t),9£{t)) is bounded on [to, +00). It follows from equation 
(7.12) that 9e is globally Lipschitz continuous on [to, +00) with some constant L > 0. 
Let us prove that 

lim 9e{t) =0. (7.39) 

If not there exist a sequence tk — > +00 and a constant c > such that 9e{tk) > c for every 
k . It follows from the Lipschitz continuity of 9^ that 

m > \ 

on the intervals {tk — ■^,tk + ■ This obviously contradicts the assumption that 9 < +00 
and concludes the proof of (7.39). 

Let us consider a sequence Tk +00 such that 9e{Tk) 0. It follows from (7.12) that 
Fg{9e{Tk),9e{Tk)) — > 0. Passing to the limit we deduce Fs{9,Q) = 0, which is impossible 
since Fs{9,Qi) = Ao{9) > 0. This shows that 9 = +oo and concludes the proof of the 
lemma. □ 
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7.3. Fast dynamics. In the rest of the section we shall prove that the function 9{t) given 
in Theorem 7.1, which determines the internal variable in the quasistatic evolution, has at 

most a jump time t (E [ta,+oo). We shall show that the sequence Oe{t) converges to 9{t) 
uniformly on the compact sets of [to,T) U (r, +oo) and that on each of the two intervals 
[to,T) and (r, +oo) 6{t) is a solution of the equation of the slow dynamics 

Bo{e{t))9{t) = Ao{e{t)) . 

Moreover Bo{d{t)) > for every t G {to,T) U (r, +oo). More precisely we shall prove that 
if Bo{9) > for every 9 G [6*0, +oo) , then there is no jump point and therefore 9{t) is the 
unique solution to the Cauchy problem 

B„{9{t))9{t)^A„{9{t)), 9{to)^eo. 

We shall prove (see Lemma 7.4) that for some values of /i we have Bo{6) > for every 
9 > 0. For other values of n the function Bq has exactly two zeros a and (3 and Bo{9) > 
for 9 e (0, a) , Bo{9) < for e {a, (3) , and Ba{9) > for e +oo) . In this case there 
is a jump time r only when 9q < j3 . More precisely if G [o^j /?) j then r = to , while t > Iq 
when ^0 S (0, a) . In both cases the jump of 9 can be determined from the relation 

9{r+) = ^{9{T-)), 

where the transition function $ is defined, through the equation of the fast dynamics 

w'{9) = B{9,w{9)). (7.40) 

More precisely for every 7 G [a, (3) $(7) is defined as the first zero greater than 7 of the 
solution to the Cauchy problem 

w'^{9) = B{9,w^{9)), w^{-i)=Q. (7.41) 

Note that 9{t—) = 9q if 9q e [a,/3) and hence t = to, while 9{t—) = a if ^0 € (0,Q:). 
The following lemma is concerned with the changes in sign of the function Bq . 

Lemma 7.4. Let 

79\/T9-344 I ^ ^/2 
/io :^ V7 + 2^19 = 0.0129... and ao := "Vv^ " ^ = 0-8639 . . 

The following properties hold: 

(a) if jjL > j2o , then Bq{9) > for every 9 > 0; 

(b) if iJb = jjbo, then Bq{6) > for every 9 ^ ao; 

(c) if jJL < jJL{) , then there exist a and (3, with 0<Q!<q;o</3< +00, such that 
Bo{9) > for every 9 G (0, a) U (/3, +00) and Bo{9) < for every 9 G {a,/3); 
moreover we have -Bq(q:) < and -Bq(/3) > 0. 

Proof. The results follow from the fact that the function 

V'{9)^V"{9) _ 1 9^ 



V'{9f-1 2(4 + 36i2)(l + 6i2)3/2 

has a positive derivative in (0,q;o), negative derivative in (ao,+oo), tends to for 6^0 
and 9 +00 , and takes the value 2yLto 3± 9 — ao- □ 

The following lemma provides some properties of the solutions of the equation of the fast 
dynamics (7.40). 

Lemma 7.5. Assume 11 < /io . For every 7 G [a,f3) there exists a point $(7) > /? such that 
the solution w-^ to the Cauchy problem (7.41) satisfies Wj(9) > for every 9 G (7,5^(7)) , 
Wj{^{'^)) = 0, and w'^{^{'-f j) < 0. Moreover the function [a, (3) (/3, $(a)] is continu- 
ous, decreasing, and 3>(7) ^ (3 as j ^ /3 . 
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Proof. By direct computations we find that 

B{0, -V'{0)) = -V'\0) > , (7.42) 

which impHes that w(^) := —V'{9) is a solution of the equation. By uniqueness we have 

that w^i9) < V'ie). 

For t; = we have that 

B{9,0) = -Bo{e)>0 for0e(a,/3), 

B{9,0) = -Bo{9)<0 for e e (/?, +oo) . ^' 

Since B{9,v) is a polynomial of the third degree in v , the sign of its coefficients implies 
that for every ^ > there exists a unique point v^{9) such that v i-^ B{9,v) is increasing 
on [O,t;«(0)] and decreasing on +oo) . It follows from (7.43) and (7.42) that 

B{9, v)>0 for 61 G (a, /3) and v e [0, -¥'(9)] . (7.44) 

Using the definitions (7.18)-(7.22), we can prove that 

lim B{9,v) = ij{-^+v + 6v'^ -Av^) (7.45) 

9— ►+00 

lim dyB{9, v) = ^(1 + 12?; - 12?;^) (7.46) 

uniformly for every v £ [0, 5] . It follows that there exists 9 > P such that 

d„B{9,v)>y (7.47) 

for every 9 > 9 and every v £ [0, ^] . 

For every v G (0, —V'{a)) let 'w{9, v) be the solution of the Cauchy problem 

w' {9,v) = B{9,w{9,v)) , w{a,v)=v, 

where prime denotes derivative with respect to 9. We want to prove that for every v G 

(0, —V'{a)) there exists 9 > a such that w{9,v) = 0. 

If not, there exists vq G {0,—V'{a)) such that w{9,vo) > for every 9 > a. By 
comparison we have also that < w{9,v) < —V'{9) for every 6 > a and every v G 
[vq, —V (a)) . Note that this inequality ensures that 9 w{6,v) is defined on [a, +00). 
Taking the partial derivative with respect to w , we obtain 

{dyw)'{9, v) = dyB{9, w{9, v))dyw{9, v) , dyw{a, v) = 1 . 

By (7.47) we have 

dvw{9, v) > dyw{9, t;)e^''(^-^") 
for 9 >9 . Integrating this inequality with respect to v between vq and —V'{a) , we obtain 

-V'{9) - w{9, Wo) = w{9, -V'{a)) - w{9, vq) > 
> {w{9, -Via)) - w{9, wo))e^'^(«-«) = {-V'{9) - w{9, wo))e^''(^-^") . 

As —V'{9) — w{9,vq) > by uniqueness and —V'{9) 1/2, we conclude that w{9,vo) — * 
—00 as 9 ^ +00, which contradicts our assumption. 

By (7.44) we have w^{9) > for 9 G (7,/?] whenever 7 G (a,/3). The same conclusion 
can be obtained also when 7 = a, using the fact that w'^{a) = B{a,0) = —Bo{a) = 
and w'^{a) = deB{a,0) = —B'Q{a) > (this can be easily deduced from equation (7.41) 
and from Lemma 7.4). On the other hand by comparison we have w^{9) < w{9,v) for 
every v G (0,-^(0;)). Since the function w{-,v) vanishes at some point we conclude 
that w-y vanishes too. Then we define 3>(7) as the smallest ^ > 7 for which 'Wj{9) = 0. 
From the previous discussion it follows that $(7) > /?, hence w;'^($(7)) < by (7.43). 
This fact, together with the continuous dependence on the initial data, implies that $ is 
continuous. It follows from a comparison argument that a < 9\ < 92 < (i implies that 
/3 < $(^2) < '^{9i) < ct- Arguing as before, we can prove that the solution tiifj of the 
Cauchy problem (7.41) for 7 = /3 satisfies w'ij{(3) = and w'^{P) < 0, hence Wjs{9) < for 
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in a right neighbourhood of j3 . By the continuous dependence on the data this imphes 
that $(7) ^ /3 as 7 ^ /3. □ 

7.4. Complete description of the quasistatic evolution. We are now ready to prove 
the main result of this section. Let Ho, a, and /3 be as in Lemma 7.4. 

Theorem 7.6. Under the hypotheses of Theorem 7.1 there exists a unique approximahle 
quasistatic evolution with boundary datum w and initial condition {uq, eo,po, zq) . It is 
spatially regular and satisfies (7.2) and (7.3). 

(a) If li> fJo or if II < Ho and 9o > 3 , then in the interval (to, +00) the function d{t) 
coincides with the unique absolutely continuous solution of the Cauchy problem 

Bo{9{t))e{t)=Ao{9{t)), 0{to)=Oo. (7.48) 

(b) If fi < Ho and a < 0o < P , then in the interval {to, +00) the function 0{t) coincides 
with the unique absolutely continuous solution of the Cauchy problem 

Bo{6{t)) m = Ao{9{t)) , e{to) = ^9o) . (7.49) 

(c) if fi < Ho and < 9o < a, then there exists a unique t > to with 0{t) = a, 
such that on the interval {to,T] 9{t) coincides with the unique absolutely continuous 
solution of the C'a,uchy problem 

Bo{9{t))9{t)^Ao{9{i)), 9{to)=9o, (7.50) 

while on the interval (r, +00) the function 6{t) coincides with the unique absolutely 
continuous solution of the Cauchy problem 

Bo{6{t)) 9{t) = Ao{9{t)) , e{T) = $(a) . (7.51) 

RemEirk 7.7. In all cases 



V3 1 



e{t) +00 and t - ip{t) 
as t +00 . By (7.2) we have 

as t +00. As 9{t) is increasing we deduce from (7.2) and (7.3) that \(T{t)\ is decreasing 
for t > to, reflecting the fact that the material softens as plastic deformation proceeds. The 
asymptotic value of the stress coincides with the value that can be obtained in the perfectly 
plastic case with elastic domain 

i^oo := n K{-V'{0)) , 

where K{Q is defined by (2.14). This reflects the fact that, with this loading program, 
the material behaves in the weakest way permitted by its internal variable as time tends to 
infinity. 

In case (a) the functions e, p, and z are continuous with respect to t . In case (b) they 
are discontinuous only at t = to while in case (c) they are discontinuous only at a point 
T > to (see Figures 1, 2, 3). 

In case (a) the energy inequality (5.13) of Remark 5.5 holds with equality. In cases (b) 
and (c) the inequality is strict when T is larger than the discontinuity time. In other words, 
an instantaneous dissipation occurs at the discontinuity time. 

Proof of Theorem 7.6. We begin by studying case (c). Let (f: (0,a) U (/3, +00) (0, +00) 
be the C°° function defined by 

(T.52) 
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Figure 2. Stress vs. imposed strain in case (b). 
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Figure 3. Stress vs. imposed strain in case (c). 



By direct computations we see that 




lim ip'{0) = 0. 
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For every e > it is useful to introduce the £ -scaled version of the velocity, considered 
as a function of 9, and denoted by Vs{0) . It is characterized by the equality 

v,{ee{t)) ■= eeS) , (7.56) 

and satisfies the condition 

where prime denotes the derivative with respect to 9 . From (7.12)-(7.22) it follows that 

v'M=e^ + B,{0,v,{0)), (7.58) 

where 

B^{9,v) := B{e,v) - eAi{e,v) + eA2{0,v)v . (7.59) 
By (7.34) and (7.38) wc have 

<VeiO) < -V'{e) for every 6* > 6*0 . (7.60) 

The proof will be split into several steps which describe the behaviour of 6e{t) in different 
intervals. Since we are considering case (c), we have fJ. < Ho and < 9o < a. 

Step (c)i . Behaviour near Oq. Let us fix 6i > 6*0 such that Bq{6) > for 9 E [6*0,6*1]. 
Let < ?7i < min{(^0)^i — ^0,1} 1 and let ei e (0,1). For every 77 G (0,r/i) and every 
£ e (0, £1) let Oj .^ be the the largest element of (^0,^1] such that Vs{6) < s{f{0) — r/) for 

every 0e [00, 

In the interval {9o,0\^^] it is convenient to write equation (7.58) in the form 

v'M = - Bo{e) + R,{e, v,{e)) , (7.61) 

Ve[9) 

where, by (7.18) and (7.59), 

R,{e,v) :=-eAi{e,v)+sA2{e,v)v + Bi{e)v + B2{e)v'^-B3{e)v^. (7.62) 

By (7.24)-(7.29) and (7.60) there exists a constant ro, independent of £ and r], such that 
the rest Rc{9,Vs{9)) can be estimated by 

\Re{9,v,{9))\<ros (7.63) 

for aU 9 G [6*0, 6*1^]. Since 

^ > 1 



v,{9) - ip{9) - r, 
for every 9 e (^o,^^,,,]) recalling (7.52) we obtain 

vm>^^^B^{9) + R,{9,vM) 

in the same interval. By (7.54) and (7.63) for every 9 G (^o,^e,r,] we have 

Vg(9) > — ri - ros, 

where 6o > is the minimum of Bq{9) on the interval [^c^i]. This implies 

Ve{9)>^^{9-9,). 

whenever < £ < bot] / {2roipo) . 

Let ipi be the maximum of (p{9) on the interval [9o,9i]. As Ve{9l^^) < s(fi{9l ,^) < e(pi, 
we obtain 

9l„-9o<'-^-^<9,-9o, (7.64) 
' ' V bo 
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whenever < £ < ei(r?) := min{?75o/(2roV3o), ?7((?i — 9o)bQ/{2(po(pi),si} . This impHes that 

v,{elv) = <viOiv) - V) ■ (7.65) 

Let tl^^ be the time such that 0siti^) = 91^. From (7.23), (7.26), (7.61), and (7.63) we 
obtain 

^^ao <<(^) + 2/x + roe in(^o,^y, 
where ao := -v^l^gl^' (6*0) > 0. This imphes 

tin -to = e f " < 1 ^;;(0) dO + - 0o) + "Mel,, - eo) < 

^ ^ Cdi \ _L ^ ^l^'Pofi . £^ ^ro<fo<^i 

S fy^e n) "I T 1 T • 

oo '' rj ctoOo V aobo 
Hence there exists a constant ci such that 

tlr,-to<ci^ (7.66) 

for < < r/i and < e < £1 (77) . 

Step (c)2 . Behaviour between 60 and a. We want to prove that for every rj e (0,r/i) there 
exists £2(??) S (0,£i(77)) such that 



-m <V (7.67) 

for every s G (0,£2(r?)) and every 9 G [61,^, a — r]]. Since this inequality is satisfied for 
9 = 91,^ by (7.65), considering the equation satisfied by Vg/s, it is enough to prove that 

Bo{9) + Re{9, e{<fi{e) + r,)) < e^'{0) < -4§r^ ' Bo{0) + Re{e, £(</^(^) - v)) ■ 



Arguing as in the previous step, we can show that Re{9,e{ip{9) ± ry)) — » as £ — » 0, 
uniformly on [9o,a — rj\. The conclusion follows from the fact that eip'{9) uniformly 
on the same interval and that 

Boie) = -^^Bo{9) < , -4§^ - Bo{0) = -j^BoiO) > 



on [6*0, a — ??] . This concludes the proof of (7.67). 

Let t1^^ be the time such that 9e{t1,,^) ~ a — rj. By (7.67) we have 

\9S) - ^mm < V (7.68) 

for every t G [i^,,,, i^,,,] • From (7.55) it follows that for every t G [ig,,,, ig,,,] 

< V^o - V < Osit) < ip]!> (7.69) 
where ipo is the maximum of (p{9) + 77 on the interval [9o,a — rj]. This implies that 

a-V- 0l,r, 2 _ .1 ^ ^ - ^7 - ejr, 
— ''s.n ''s.n — 

<Po-V 

By (7.66) for every t] G (0,r/i) we have — > to as £ — > 0. By compactness we can 

assume that there exist t"^ > to such that t^ ^ ^ t^-^ as £ ^ along a suitable sequence. 
As tl^^ < t"^ g foY S < T] , we have t"^ < tg for S < i]. This implies that there exists 
r e (to, +00) such that t^ ^ r as 77 — > 0. By (7.69), using a diagonal argument, we can 
find a subsequence, still denoted 9e, and a locally Lipschitz function 6^: (to,r) — > M. such 
that for every < 77 < Tji and every to < < ^2 < t"^ we have 

9^ -> 9^ uniformly in [ii,i2] , (7.70) 
9e L weakly* in L~([ti, £2]) . (7.71) 
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Passing to the limit in (7.68), we deduce that — (p(6'*(t))| < rj almost everywhere in 

{to, t'^)- By the arbitrariness of > we get that 6^ is a solution of the equation 

e^t) = ^{e,{t)) (7.72) 

for t e {to,T). 

Let us fix e (0,?7i). By (7.66) for every t G {to,t^) we have that t G (i^,,,, *£,,,) for £ 
small enough, which gives, thanks to (7.69), 

|^e(0-C,l<(*-^e,>S<(^-*0)^3- 

Passing to the limit as e ^ 0, and using (7.64), we deduce that 

\0*{t)-9o\ < (t-io)<Po- 

This shows that 9^, satisfies 9^,{t()) = 9() . By uniqueness of the solution of the Cauchy 
problem (7.50) wc deduce that 9^,{t) — 9(t) for every t G {t[j,T). This implies that the full 
sequence 9^ tends to as e ^ 0, uniformly on compact subsets of {tQ,T). 

As & {to,T) and 9e{t1^^) = a — r], we deduce that 9{t'^) = a — r], and hence 

9{T)=a. 

Step (c)3 . Behaviour near a. Assume rj € (0, ?7i) and < e < miD.{r]/ {(p(a — r])+ri), e2{ri)} . 
Remark that (7.67) implies Ve{a — r)) <r]. We want to find 773 e (0, rji) such that for every 
ri e (0,%) there exists £3(77) € (0, £2('?)) for which 

Ve{9) < V (7.73) 

for every e G (0, £3(77)) and every 9 G [a — T],a]. Let 9 bo the largest element of [a — rj, a] 
such that (7.73) is true for every 9 G [a — r],9]. By (7.62) the rest of equation (7.61) can be 
estimated for e < rj < rji and 9 G [a ~ rj^O] by 

\Re{9,v,{9))\<nr, (7.74) 

for a suitable constant ri independent of 77 and e. Since Bo{9) > for 9 G [a — r],9], by 
(7.23) and (7.61) we have 

where clq := n\£,o\. 

Let Ws be the solution of the Cauchy problem 

= 73Y w{a - T]) = aoV^ ■ 

By (7.67) for e small enough we have Vs{a — rf) < e{ip{a — ri) +ri) < CbQ^/e, which implies, 
by a comparison argument, that 

Ve{9) < w,{9) (7.75) 

for every 9 G [a — ri,9\. Since is increasing, we have w'^{9) < yje + rirj, which gives 

We{9) < do^/e + {^/e + rir]){9 - a + rj) . 

If T) < 7^3 := min{l/ri,77i} and e < £3(77) := min{(l - 7?i)^7?^/(ao + ri)'^ , r] / {ip{a - rj) + -q), 
£2(77)}, by (7.75) we have Vs{9) < 77 for every 9 G [a — 77,^], which implies that 6 = a and 
concludes the proof of (7.73). 

Let tl be the time such that 9e{tl) = a. Let us prove that — t^^^ is uniformly small. 
To this aim, we observe that (7.23), (7.61), and (7.74) give 

v'e{9) > - Bo{9) - ri?7 for every 9 G {a - r],a) , 

Ve[9) 
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where qq ■= —V5ii\^q\V' (Oq) > 0. Using (7.26), (7.73), and the inequahty rirj < 1, we 
obtain for < r/ < rys and < e < £3(7/) 

^'-^Ir, = r 7r(mde<^ r U{e)+Bn{e)+rr,i) de < ^^r?. (7.76) 

Let us prove that 



Hmt =T. (7.77) 

£—►0 

For < 77 < 773 and < £ < £3(77) we have 

\tl-T\ < \tl-tlj + \tl^-tl\+\tl-T\. 

Using (7.76), we deduce that 

,3 ,211 + 2 ,2 , 

limsup \tl -t\< — 77 + \tl - t\ . 

e^o ao 

Passing to the limit as 77 — > 0, we obtain (7.77). 

Step (0)4. Behaviour between a and $(0;). Let Wa be the solution of the Cauchy problem 
(7.41) with 7 = a. We want to find 774 € (0,773) such that for every 77 € (0,774) there exist 
£4(77) G (0,£3(7/)) for which 

Wa{0) - 77 < v,{e) < Wc^iO) + 77 (7.78) 

for every s G (0,64(77)) and every 9 G [a, $(a) — 77]. This will be done by comparing the 
equations (7.41) and (7.58) satisfied by Wa and Vg, respectively. 
By (7.24), (7.25), and (7.59) we have 

B{0, v) - ais < B^{e, v) < B{e, v) + a^e (7.79) 

for 61 > 00 and t; e [Q,-V'{e)\, were ai := 4/x|^g| and 03 := -/;i|^^|/y'(6'o) . For every 
p e ]R let Zp be the solution of the the Cauchy problem 

z'p{6)=p + B{e,Zp{e)), Zp{a)=p (7.80) 

By comparing (7.58) with this equation, we obtain that, if < £ < —p/a\, then Vs{0) > 
Zp{9) for 9 > a. By the continuous dependence on initial data Zp Wa uniformly on 
[a, $(a)] as p ^ 0. Therefore, for every 77 e (0,773) there exists ^1(77) > such that 
Zp{9) > Wa{9) — rj for every p e [— pi(7/),0] and every 9 e [a, $(0;)]. It follows that, if 
77 G (0, 773) and e < £4(77) := p\{ri)/a\ , then the first inequality in (7.78) is satisfied for every 
9 e [q, *(«)]. 

In the previous step we have proved that Vg (a) < 77 for 7/ < 773 and £ < £3 (77) . Let us 
fix 77 e (0,773), p e (0, 1), and e £ {0,63(7])) with e < pi]. If v^{a) > e/p, we set 9^^ := a. 
If Ve{a) < e/p, let 9'^ p be the greatest element of [a, (3] such that Ve{9) < e/p for every 
e e [a,9lp]. Prom (7.44) and (7.79) it follows that 3^(6, v) > -aie for 9 e [a,/3] and 
V e [0,~vi9)]. Using (7.23), (7.58), and (7.60) we deduce that v'^{9) > e{an/v^{9)) - a^e , 
with ao := -\/3|Col V^'(6'o) > 0, hence Ve{9)v'^{9) > a^e - axe^ / p for every 9 G [a, 6*^,^] . 
If £ < e\{r\,p) := min{pao/(2ai), /977, £4(7;)} we obtain Vg{9)v'^{9) > eao/2 for every 6 G 
[a, 9^ p] , hence, by integrating, 

Ve{9)'^ > aoe{9 - a) for every 9 G [a, 9lp] . (7.81) 

It follows that 

0<9tp-a<^. (7.82) 

In particular, if Vs{a) < e/p and e < £4(77, p) := min{(/3— a)aop^, £^(77, p)} , we have 9f p < (3 
and 

ve{9ip) = (7.83) 
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and 

Ve{0) < V for every 6 e [a, O^p] . (7.84) 
Let be the time for which e^ir^p) = 0^p. Using (7.81) and (7.82) we deduce that 

tl = £ < =-^./6'4 -a< <— n. (7.85) 

Suppose that ij < tjI := min{ao/(3ai), «;q,(/?)/3, 773} , p G (0,1), and e < e'^{r],p) := 
mm{pril, e^{ri, p)} . Let be the greatest element of [a,P] such that 

v,{e) < 2rjl for every 9 G [a, 9*] . (7.86) 

By (7.84) we have 9'^^p < 9"^ . By the first inequality in (7.78) we have Ve{P) > Wa{l3) - t] > 
2r]l , hence 6^ < (i and 

vM) = H . (7.87) 

We want to prove that 9"^ — a is larger than a positive constant, independent of e. Using 
(7.23), (7.58), (7.79), and (7.60) we deduce that v'^{9) > s{ao/ve{9)) - eai > £((ao/(2r/l)) - 
ai) > for every 6 G [a,6f]. This implies that is increasing on [a,^^], therefore 

ve{e) > - for every 9 e [9^ 9^] . (7.88) 
p 

As Vs{a) < 7] < rjl and Vs{9f) — 2r]l, there exists a unique point 9j in {a, 9"^) such that 
v^{9^) = r]l. It follows that Vei9) > r/l for every 9 € [9t,9i]. By (7.23)-(7.29) there 
exists a constant 61 > such that B^{9,v) < bi for every e G (0,1), 9 G [60, +00), and 
V e [0,-V'{9)]. Let ao := iJ,\Q\. Using (7.23), (7.58), and (7.60) we obtain 

v'M<^ + bi<s^+h<2h 

Ve{9) Vl 

for every 9 G [^e,^^] and every e < e\{ri,p) := min{?7]&i/ao, £|(?7, p)} . Therefore (7.87) 
gives 



which implies 



9t-a>9l-9l>^. (7.89) 
2oi 

For every rj < tjI := min{/3 — a, $(a) — /3, 774/ (26i), rjl} we define 

TUri '■= Tahi{wa{9) : a + 7] < 9 < ^(a) — 77} > . 

Since ^0 as rj ^ 0, there exists 774 G (0,7;|) such that rrirj < 773 < 1 for every r] < r^. 
Using (7.89) and the first inequality in (7.78) (with 77 replaced by mjj/2) we obtain that 

Ve{9)>'^>- fov9G[9tMa)-Tl] 
I p 

provided that 77 < 774, p < 1, and e < £4(7;, p) := min{£4(77T.,,/2), p77i,,/2, £4(77, p)} . Using 
also (7.88) we obtain 

Ve{9) > - for every 9 G [9^^, $(a) - 77] . (7.90) 
P 

Let t^ jj be the time for which 9^{t'^_,^) = $(a) — 77. From the previous estimate we have 
- = e — ^ < p($(a) - a) < r7($(a) - a) , (7.91) 
if we have also p < rj. 
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By (7.23), (7.58), (7.79), (7.60), and (7.90) wc obtain 

v'M < {ao + a2)p + B{e,v,{e)) for every e [(?^^^, $(a) - 77] , 

whenever 0<?7<r74, < p <ri, and < e < £4(77,^) (recall that £4(77, p) < p). Replacing 
p by p/{ao + 02) we obtain 

v',{0) <P + Bi0, v.m for every 9 € [0^^, $(a) - r,] , 

whenever 0<r7<?74, 0<p< P2iv) ■= minlry, 77(00 + 02)}, and < £ < el{ri,p) := 
min{£|(77, p), £^(77, p/ (oq + 02))} . 

If < £ < £t(77,p) := min{£3(p,p),£|(,7,p)}, by (7.84) we have v,{elp) < p. By (7.44) 
the solution Zp of (7.80) is increasing on [a, f3] , hence Zp{9^ p) > p. By comparison we have 

Ve{d) < Zp{9) for every 9 e [9^^, $(a) - 7/] . (7.92) 

On the other hand by the continuous dependence on the data Zp Wa uniformly on 
[a, $(a)] as p — > 0. Therefore, for every 77 < 774 there exists P4{ri) G (0,^3(77)) such 
that Zp{9) < Wa{9) + T] for every p e {0, psir])) and every 9 e [a,^{a)]. By (7.92) 
we have Vg{9) < Wa{9) + rj for every 9 G [^e ^4(17)1 ^C*^) ^ v\ when rj < r]4 and e < 
£4(77) := min{£4(77, ^4(7/)), 77^} , which proves the second inequaUty in (7.78) on the interval 
p^(^), ^(a) — v] ■ The same upper bound on the interval [a, 6'^ p^(ri)] follows from (7.84). 

Step (c)5 . Behaviour near $(a) . By (7.78) we obtain in particular 

Wa($(a) - 77^) - 77^ < We($(a) - 77^) < Wa($(a) - 77^) + ?7^ 

for every 77 e (0,774) and every e £ (0, £4(r7')) . As w'^{^{a)) = —Bo{^{a)) < 0, there exist 
two constants 02,03 > and a constant 775 e (0,774) such that 

C2V'^ < We($(a) - 7?2) < C37?2 (7.93) 

for 77 e (0,775). Note that for e < £5(77) := min{c277^/(iy?($(ci!) - 77^) + 77), £4(77^), £4(77^)} we 
have 

e{ip{^{a) ~ 7f ) + ri) < Ve{^{a) - jf) . 

For any r] < ri\ and £ < £5(77) let be the largest element of [$(0;) — 77^,^(0!) + 1] 
such that 

e{if{9)+r,)<v,{9)<czri' 

for all 9 G [i>(a) - ri^,9l^^]. Using (7.16)-(7.29) the rest Re in equation (7.61) can be 
estimated by 

\Re{e,Vem\<r2{e + ri), (7.94) 
for every 6 G [^(a) — where r2 is a constant independent of 77 and s. Since 

£ 1 

MO) - f{9)+v 

for every 9 € [^{a) — V^, 9^^^] , recalling (7.52) we obtain from (7.61) 

v'A9) < -^l^B,{9) + Re{9,vM) 

in the same interval. Let 60 > be the infimum of Bq{9) for 9 > ($(a) + /3)/2, and let 
(po be the supremum of (p{9) + 1 on the same half-line. By (7.54) and (7.94) in the interval 
[$(q) — 77^, 9^^^] we have the estimate 
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provided that r] < rj^ := min{6o/(4r2<^o), ??5} and < £ < £| := inin{77, £5(77)} . As 
Ve{0) > by Lemma 7.2 and by (7.56), we obtain 

where the last inequaUty follows from (7.93). Therefore there exists % € (0, Tjf) and for 
every -q S (0,775) there exists £5(77) G (0,£4(7/)) such that 

< - *(a) + 772 < ^$^77 < 1 (7.95) 

for every e G (0,£5(77)). As v'^{6) < on [$(a) — 77^, ^e,^], from the maximality of 0^^^ we 
deduce that 

vM,r,)=<'PiOlr,) + V). (7.96) 

Let t^,^ be the time for which Oeit^,^) = 6^^^. Since Vs{6) > e{ip{6) + rj) for every 
6 e [^{a) - i!f,Ol^^] , using (7.54) we obtain that 

ft, - '1,. = . /•■■■ < /•■■■ ^ < - *(«) + V') 

which, by (7.95), implies 



<I>(a) — r;^ ^£ 

(p{e) +77 ipo 



tl,-tl,^<=r^V (7.97) 
ooVo 

for every rj € (0,7^5) and for every e G (0,£5(7/)). 

Step {c)q. Behaviour between $(a) and +00. Using (7.96) and arguing as in Step (c)2 we 
can prove that for every t] G (0,775) there exists £6(77) G (0,£5(77)) such that 

Ve{0) 



e 

g5 



< V 



for every e G (0,£6(??)) and every 6 G (f?5,p+oo). 

From (7.95) it follows that, passing to a subsequence, we have 

(^e,ri ^ as £ ^ and 9^ i>(a) as 77 ^ . 

Using (7.85), (7.91), and (7.97) we get 

< 1^2 ^ _ ^)^2 ^ 2C3^^ 

for every 77 G (0,775) and for every e G (0,£5(7;)) (recall that this implies e < £4(77^), hence 
£ < £4(77^,^3(77^)) and £ < £4(77'^, P3 (77^))). Passing to a subsequence we obtain 

^e,)7 ^ as £ — > and — > r as 77 — > . 

Arguing as in Step (c)2 , we deduce that 9e{t) converges to the solution 9{t) to the Cauchy 
problem (7.51) uniformly on compact subsets of (r, +00). This concludes the proof of 
case (c). 

In case (b) the proof is divided into three steps. 

Step (b)i . Behaviour between 9q and $(^0) • It is enough to repeat the proof of Step (0)4 
with a replaced by 6*0 ■ 

Step (b)2 . Behaviour near $(^0) • It is enough to repeat the proof of Step (0)5 with a 
replaced by 9q . 

Step (b)3 . Behaviour between $(^0) and +00. It is enough to repeat the proof of Step (c)6 
with a replaced by . 
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In case (a) we first assume that Bo{0) > ior 9 > 9o- This happens when ij,> jiQ, ov 
when n = jjLQ and 6q > ao, or when /z < /xo and 6o > (3. In all these cases the proof is 
divided into two steps. 

Step (a)i . Behaviour near Oq ■ This is identical to Step (c)i . 

Step (a)2. Behaviour between Oq and +oo. It is enough to repeat the proof of Step (c)6 
with $(q!) replaced by ■ 

In the case ji < jiQ and Oq = (3 , ot ji = jiQ and = ckq , we have only to modify the 
proof of Step (a)i . First of all we find r/i e (0, 1) such that for every r] e (0, r]i) there exists 
£i(r?) S (0,ry) for which 

Ve{0) < V 

for every e G (0,ei(77)) and every 9 G [P,f3 + rj] (with (3 — aQ for fi = iiq). This can be 
proved as in Step (c)3 , replacing a — rj by /?. 

If {P + rj) > e{(f{(3 + rj) + r]), we define 9* ,^ as the largest element of [f3 + r],(3+l] such 
that 

eMe)+v)<v,{9)<r] 
for all 9 e [(3 + ri, 9'*^^] . As in Step (c).-, , wc can prove that 9* ^^ /3 + ri a.s e ^ . 

If Vs {(3 + r]) < £{<f {(3 + ri) — r]), we define 9*^^ as the largest element of [P + ri,(3+ 1] such 
that 

Ve{e) < S{^{9) - 7?) 

for all 9 G [(3 + ri, 6** ^] . As in Step (c)i , we can prove that 0*^ ^ f3 + rj as e— >0. 
If £(¥'(/? + v) - V) < Ve{^ + v) < + r/) + r;) , we define 9*^^ := P + rj. 

This concludes Step (a)i in the case /x < /xq and 9o = (3. 

It remains to study the case jjL = fiQ and 6*0 < ao . The proof is split into four steps. 

Step (a')i . Behaviour near 6o ■ This is identical to Step (c)i . 

Step (a')2. Behaviour between 9q and a^. This is identical to Step (c)2. 

Step (a')3 . Behaviour near ao . The behaviour in the interval [ao — r], ao] can be studied 
as in Step (c)3 , while the behaviour in a right neighbourhood of ao can be studied as in 

Step (a) 1 for the case /i < /iq and 9o = j3 . 

Step (a')4 . Behaviour between ao and +oo . It is enough to repeat the proof of Step (c)6 
with p replaced by ao . □ 

8. Examples with concentrations and oscillations 

In this section we consider solutions of the reduced problem in dimension d = 1 considered 
in Section 6 and corresponding to simple shears. To simplify the exposition we take 

/i=i and := {(a,C) e MxM : |a| + Id < 2}, (8.1) 

so that 

a''{t)^e^it) and ff«(e, 0) = 2(|d V |e|) . (8.2) 
About the softening potential V we assume that 

V{9) = V{-9) for every 9 gR (8.3) 

and 

V'{e) = -1 if and only if 61 > 1 . (8.4) 

Together with the general assumptions (2.31) and (2.32), these conditions imply that —1 < 
V'{9) < 1 for every 61 G M and 

-1 < V'{9) < for < 6* < 1 . (8.5) 
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We will consider the boundary datum 

w'^{t,y):=ty (8.6) 
and the initial condition {uq, Cq^Pq, Zq) , where 



R ^ n „R - n „R - n ^„ r i i 



g — VJ , Cq — U , pQ — U Wll [ 2 ' 2 J ' 

_ _l ^ IB ;o Q 

2' 2 



p^ = on[-i,i], (8.7) 



and : [—h, h] ^ R is a, suitable piecewise continuous function. 



8.1. Strain localization. The following theorem describes an example where the strain 
has a singular component supported by a lower dimensional manifold. 

Theorem 8.1. Besides (8.1)-(8.7), assume that Zq is continuous and Zq{0) = 1 > ZQ{y) > 
for every y ^ 0. Let {u^,e^,fj,^) be a solution of the reduced problem for the approx- 
imahle quasistatic evolution with boundary datum and initial condition {uq, CqjPq, Zq) . 
Then for every t G [0, +oo) we have 



ty if tG [0,1], yG [-11, 

u«(i,y) 



y + t-1 «/ t e (1, +oo) , y e (0, i] , (8.8) 
^y-t+1 iftG{l,+oo), y e [-i,0), 

e«(i) =tAl, /xf = Vw,^«(t)) ' (8-9) 

where 

p«(i):=(i-l)+<5o, z''it):=z^ + {t-l)+6o, (8.10) 
5q being the unit Dirac mass concentrated at . Moreover, for every T > we have 

VHn{ii"-0,T) = 2{T-l)+ , (8.11) 
and the energy inequality (ev4) of Theorem 5.4 holds with equality. 

Remark 8.2. In this example we see the phenomenon of strain localization at y = even 
if the boundary and initial data are smooth. If the reduced one-dimensional problem is used 
to describe the evolution of simple shears, this example exhibits the formation of a shear 
band on the plane Xi = for every time t > 1. 

Proof of Theorem 8.1. By Remark 4.10 for every £ > the solution {u^, ef-,pf , z^) of the 
reduced e -regularized evolution problem with boundary datum and initial condition 
^Q^Po^^o) satisfies (4.76) and (4.77). Moreover zf{t,y) is continuous with respect to 
(t, y) . We observe that 

N^Knia, C) = ^ ((a + C - 2)+, (a + C - 2)+) (8.12) 

for every £ > and every (a, C,) G RxR with a > and — 2 < ^ — a < 2. 
Let te be the largest element of [0, +oo] such that 

-2<-V'{zf{t,y))-cT^{t)<2 (8.13) 

for every t G [Q,te) and every yG[-\,\\. By (4.76), (4.77), and (8.12) we have 

e^{t) = l-^J\{ef{t)-V'{zf{t,y))-2)+dy, (8.14) 

zf{t, y) = l(ef (t) - V\z^{t, y)) - 2)+ for every y g[-\, \] (8.15) 

for every t G [0, te) . 

For t G [0, 1] the functions 

ef{t)=t and zf{t,y)=zo{y) (8.16) 
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solve (8.14) and (8.15) and satisfy the initial conditions ef^(O) = and z^{0,y) = ZQ{y). 
Moreover they satisiy (8.13) with strict inequalities. Therefore te > 1 and (8.16) gives the 
solution of the reduced e-regularized evolution problem for t e [0, 1] . 
We observe that 

z^{t,y)>z^{y) (8.17) 
for every t € [0,t^) and every y G [— |, |] . This follows easily from (8.15) and from the 
initial condition z^{0,y) = ZQ{y). 
Let us prove that 

ef{t) > 1 (8.18) 

for every t € (l,ie). Indeed, this is true for t near 1, since ef{l) = 1 and e^(l) = 1. 
If (8.18) is not satisfied for every t G (l,ie), we can consider the first point Te > 1 where 
ef(re) = 1. It is clear that ef(re) < 0, but from (8.14) we obtain ef(re) = 1. This 
contradiction proves (8.18). 
Let us prove that 

lim sup ef(t) = 1. (8.19) 

i<t<t. 

Since the function x V'{zq{x)) is continuous and Zq{0) = 1, for every tj > there 
exists ^ > such that 

-V'{z^iy))>~V'{l)-r] = l-v, 
whenever \y\ < 6. Together with (8.17) this imphes that 

-V'iz^{t,y))>l-rj, 

for 1 <t <t£ and |a;| < 6. Therefore, (8.14) implies 

ef (i) < 1 - ^ f^iefit) - V'izf{t,y)) -2)+dy<l- ^(ef (i) - 1 - r?)+ . 

As (1) = 1, by a comparison argument we obtain 

sup e^{t) <l + r]+ — . 

l<t<ts 

This implies 

Umsup sup ef^{t) < 1 + 7] , 

S^O l<t<te 

which, together with (8.18), gives (8.19) by the arbitrariness of r]. 
We deduce from (8.16), (8.18), and (8.19) that 

ef(t)^iAl (8.20) 

uniformly with respect to t > 0. It follows from Definition 5.1 that the first equality in (8.9) 
is satisfied for every t G O , and hence for every t G [0, +oo) by left-continuity. 

As z^{y) > 0, by (8.5) and (8.17) we have -V'{z^{t,y)) > 0. Therefore, by (8.19) 
there exists eo > such that -3/2 < -V'{z^{t,y)) - af{t) < 1 for every e € (0,eo), 
t e [0, fe), and y G [—5,5]- By continuity and by (8.13) this implies that = +00 for 
every e € (0, £0) . 

For every e G (0,£o) and every t > 1 we define 

aS) ■■= snp{\y\ : zf{t,y) > z^{y)} . (8.21) 

Since z^{t,y) is increasing with respect to t by (8.15), we have that a^it) is increasing 
with respect to t. Let 

ae(+oo) := lim ae{t) . (8.22) 
t— ^+00 



Let us prove that 



limae(+oo) = 0. (8.23) 

£— >0 
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Let us fix ?7 > 0. As ^^(y) < 1 for \y\>ri, by (8.5) 

max~r{z^{y))^l~Sr, (8.24) 
\y\>n 

for some 5^ > 0. By (8.19) there exists € (0,£o) such that e^{t) < 1 + Srj for every 
e < Er, and every t>l. Therefore, if \y\ > r], by (8.24) we obtain 

{ef{t)-V'{zl^{y))-2)+=0 

for every e < £^ and every t > 1. This imphes z^{t,y) = ZQ{y) by the uniqueness of the 
sohition of (8.15) with the initial condition ZQ{y). We deduce that ae{t) < rj for every 
e < En and every t>l, which gives (8.23). 

To continue the proof of the theorem we need the following lemma about concentrations 
of positive functions or measures. 

Lemma 8.3. Let U he a hounded open set in , d > 1 , and let pk he a sequence in 
M+(U) which converges to p weakly* in Mi,{U) . Assume that p is singular with respect to 
the Lebesgue measure. Then Sp,. — ^ 5p weakly* in GY{U; R) . 

Proof. Passing to a subsequence we may assume that 6p^ /j. weakly* in GY{U ; K) . Then 
p = bar(/i) by (2.5). As suppJp^ C J/x [0, +oo)x [0, +oo) , we have also 

supp/x C C7x[0,+oo)x[0,+oo). (8.25) 

Let us consider the representation of given by [6, Theorem 4.3 and Remark 4.5]: 



(/,m)=/ f{x,^,l)dij?'{x,0+ [_ /(x,^,0)dM°°(^,O = 

Juxm. "'c/xSm 

= / ([ f{x,^,l)df,^'''{0)dx+ nf f{x,^,0)df,^'°°{^))dX°°{x) 



where ji^'^ and arc probability measures and \°° ttjj{ii°°) . By (8.25) we have 

supp/i^'^ C [0, +oo) and supp/Lt^'°° = {1}, hence = Si, the unit Dirac mass at 1. 
Let 

f+OO p 

u'-ix) := / erfM"'^(0 , u°°{x) := / ^^/."'^(O = 1 • 

Jo JEr 

By [6, Remark 6.3] we have 

p = bar(/i) =u^ + w°°A°° = + A°° . 

Since p is singular, this implies = a.c. on U and A°° = p. As /U^'^ is a probability 
measure, it follows from the definition of that H^'^ = So. Therefore, (8.26) yields 

(/,m)= / f{x,0,l)dx+ [_f{x,l,0)dp{x), 
Ju Ju 

which, by (2.6), is equivalent to fi = Sp. □ 

Proof of Theorem 8.1 (Continuation) By (4.2) and (8.12) for e < eo we have zf{t,y) = 
p^{t,y), so that the initial conditions (8.7) give 

z^{t,y)=p^{t,y) + z^{y) (8.27) 

for every t>0 and every y e [— ^, ^] . In particular (8.17) implies that 

pf(t,2/)>0 for every ?y e [-i,i], (8.28) 

while (8.21) and (8.22) yield pf{t,y) = for \y\ > a^{+oo). We deduce that for every 
t > there exist a sequence ej — > and a constant (p{t) > such that (t) ^ 'p{t)So 
weakly* in Mb([-i,i]). Since ef.{t) ^tAl = e^{t) and {uf^{t),ef:{t),pf.{t)) satisfies 
condition (cvl)^ of Definition 4.1, arguing as in the proof of Theorem 5.4 we obtain that 
u^.{t) converges weakly* in i?y([— i,i]) to a function u*{t) S BV{[—^,^]) such that 



QUASISTATIC EVOLUTION IN PLASTICITY WITH SOFTENING 53 

{u* {t),e^^{t),ip{t)do) satisfies condition (evl)^ of Theorem 5.4. Using this property and 
(8.6) we deduce that u*{t) coincides with the right-hand side of (8.8) and (p{t) = {t— l)"*" . 
As the limits do not depend on the sequence Sj , we obtain that 

u^{t) u*{t) weakly* in BV{[-^, i]) , (8.29) 
Pf{t) ^ weakly* in Mfc([-i, i]) , (8.30) 

zf(t)-z«(t) weakly* in Mb([-i,i]), (8.31) 

where p^{t) and z^{t) arc defined by (8.10). It follows from Definition 5.1 that (8.8) is 
satisfied for every t € Q, and hence for every t G [0, +oo) by left-continuity. 

Let V: [-i,i]xRxM-^ [-i,i]xMxMxM be defined by ^(?/,/3,r?) := {y, 0, P+r]z^{y),'n) . 
It follows from (8.27) that 

^(p«(t),z«(t)) = V'(^pR(t)) , ^(pR(t),2«(t)) = V'(^pH(t)) • (8-32) 

By (8.28) and (8.30) we obtain that 5pK(t) ^ dpR(^t) weakly* in GY{[-^, ^];R) thanks to 
Lemma 8.3. By (8.32) we conclude that 

'5(pH(t),2«(t)) ^ ^(p«(t),z«(t)) (8-33) 

weakly* in GY{[—^, i];Rx]R). It follows from Definition 5.1 that the second equality in 

(8.9) is satisfied for every t € 8, and hence for every t S [0, +00) by left-continuity. 
To prove (8.11) it is enough to show that that for every ti and ^2, with < ti <t2, 

(iJ«(/32 -/3l,e2 - 0l),/X(^(,(j/,/?l,^l,/32,^2,7?)) = 2{{t2 - 1)+ - ih - 1) + ) . 

Taking into account (8.2) and the left continuity, we have to prove that 

{\p2 - V 1^2 - 6i\,tif^t^{y,puOi,p2,e2,v)) = {t2 - 1)+ - (ii - 1)+ (8.34) 

for every ti, t2 & Q, with < ti < t2- We know that for every ti, t2 & @, with < ti <t2, 
there exists a sequence ej — > such that 

'^((pfl.(tl),Z«.(tl)),(p« (t2),Z« (t2))) ^ Mt1t2 ' 

hence 

(1/32 - /3i| V 102 - e,\,ti^^t^iy,pi,euP2,02,v)) = 

1 

= lim |pf (i2,y)-pf (ii,y)| V|z^.(t2,y)-z^.(ti,2/)|rf2/ = 
= .1™ /' {p^^it2,y)-p^^{h,y))dy, 

J— >00 J_l ^ J 3 ' 

where in the last equality we used (8.15) and (8.27). Therefore (8.30), together with the 
definition of p^(i) given in (8.10), yields (8.34). 

The equality in (ev4) of Theorem 5.4 follows now from (8.9), (8.10), and (8.11). □ 

8.2. Oscillation of the internal variable. The following theorem describes an example 
where strain localization occurs together with a strong oscillation of the internal variable. 

Theorem 8.4. Besides (8.1)-(8.7), assume that the Junction Zq is odd on [—5,5], con- 
tinuous on [— ^,0) U (0,^], and satisfies Zq{0+) = 1 > -2(f(y) > for every y ^ 0. Let 
(it^, e^, fj,^) be a solution of the reduced problem for the approximable quasistatic evolution 
with boundary datum and initial condition {u§, e§, p^,z^). Then for every t G [0, -|-oo) 
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have 



ty i/ie [0,1], ye 

y + t-1 t e (1, +00) , y e (0, i] , (8.35) 
y-t+1 if tG{l, +00), ye [-1,0), 



e^{t) -tAl, Mf - i<5(pH(t),z«+pH(t)) + 5^(p«(t),z«-p«(0) ' (8.36) 



where 



p^{t) := {t - l)+5o , (8.37) 
So being the unit Dirac mass concentrated at . Moreover, for every T > we have 

VHnip."; 0, T) = 2{T - 1)+ , (8.38) 
and the energy inequality (ev4) of Theorem 5.4 holds with equality. 

Remark 8.5. Let , C" : [0, +00) ^^([-i, i]) be the functions defined by z^{t) := 2^ 
and C^(0 := -V'iz^) for every t > 0. It follows from (8.36) and (8.38) that 

C'it) = -7r^({r }/xf ) and (p«(i), z^it)) = bar(Mf ) 

for every t > 0, so that in this case the stress constraint (ev3) of Theorem 5.4 depends on 
the Young measure /j,^ only through its barycentre. However, the barycentres and 
do not satisfy the energy inequality (5.13). Indeed, if T > 1, by (8.37) we have 

Ml|e''(T)||i = i , VMp", ^""l 0, T) = 2(r - 1) , 

V{z^{0)) = V(z«(T)) = V(z^) , {a^{t),Dw'^{t)) = t A 1 , 

which contradicts (5.13). Note that by (8.36), (8.37), and (8.38) we have 

VHHit^^;0,T) = 2{T-l), 

{{V}{e, v), ti§{y, b, e, rf)) = V{z^) - (T - 1) , 

so that the energy inequality (ev4) of Theorem 5.4 holds with equality. 

This example shows that, if we consider just the barycentres, which coincide with the 
weak* limits of the solutions [u^ , , , zf) of the reduced £-regularized evolution prob- 
lems, the energy inequality cannot be obtained because, by the lack of convexity, we neglect 
some important terms generated by the space oscillations of the approximate solutions, that 
are captured only by the Young measure formulation. 

As in Theorem 8.1, in Theorem 8.4 we see the phenomenon of strain localization at y = , 
even if the boundary and initial data arc smooth. In the latter we sec also, for every time 
t > 1, a strong oscillation of the inner variable z^ concentrated near the point where the 
plastic strain is localized. 

Proof of Theorem 8.4- Let Zq : [— |, |] — > M be the function defined by 

z^iy) = z^{y) for y e (0, \] , 

zo«(0) = Zo«(0+) , (8.39) 
z^(y) = -zj^{y) for y G [-\, 0) . 

It turns out that Zq is even and continuous, and satisfies Zq {Q) = 1 > z^{y) for every 
y ^ 0. Let {uf, ef,pf, zf) be the solution of the reduced e-regularized evolution problem 
with boundary datum and initial condition (uq , Cq , Pq- , Zq-) and let {uf',ef^,p^,zf ) 
be the solution of the reduced e-regularized evolution problem with boundary datum 
and initial condition {uq, eQ,pQ, Zq) . By the symmetries of the problem we can prove that 
for every t>Q and every y G [— 5, ^] we have 

uf{t,y) = uf(t,y)^-uf{t,-y), ef{t)=ef{t), 

pf{t,y)=pf{t,y)=pf{t,-y), ^ ' ^ 
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and 

z^{t, y) = -zf{t, y) = -zf{t, -y) for y e (0, \] , 
zf{t,y) = -z^{t,y) = -zf{t,-y) for e [-i, 0) . ^ " ^ 

Since Zq satisfies the hypothesis of Theorem 8.1, we deduce from (8.20), (8.29), (8.33) , 
and (8.40) that 

u^{t) u*{t) weakly* in BV{[-^, i]) , (8.42) 
ef(i)^tAl, (8.43) 
%f (0,if (t)) ^ ^ip'^W.zll+p^ii)) weakly* in GY{[-^, i]:MxM), (8.44) 

where u*{t) is defined as the right-hand side of (8.35) and p^{t) is defined in (8.37). 
Let us prove that for every t>0 

<^(pf (t),zf (t)) ^ 5'5(pH(4)^^H+pi,(4)) + ^'5(p«(t),^«-pH(t)) • (8.45) 

To this aim, we fix / G C'*°™([-i, i]xRxRxM) and observe that, by (8.40) and (8.41), we 
have 

{f,S{pf{t),z^{t))}= / J{y,pf{t,y),zf{t,y),l)dy = 

1 jL 

h{y,p^{t, y), z^{t, y), ^)dy+lj\ h{y,Pe{t. y),^^it, y), 1) dy , 



where 



My,PAv) = r'^^'''''^ '''-'^ (8.46) 
^ ' \f{-y,fi,e,v) ify<0, ^ ' 

M,,/.,M) = ({J-^/'-^^) (8.47) 

[f{y,P,-0,v) ify<o. 

Using (8.44) we get 

lijn(/,^(pf (t),zf(t))) = 5(/i'^(p«(t),f?+P«(t))) + hif2^^(pR(t),z^+pR{t))) ■ 
Since p^{t) = {t- l)+So, using (2.6), (8.46), and (8.47) we obtain 

i(/l,lJ(pH(t),J«+pH(t))) + |(/2.<5(pH(t),2«+pH(t))> = 
= 5(/'^(p«(t),Z«+P«(t))) + |(/''^(p«(t),2?-p«(t))) ' 

which concludes the proof of (8.45). 

It follows from (8.42), (8.43), (8.45), and from Definition 5.1 that (8.35) and (8.36) are 
satisfied for every t E S , and hence for every t € [0, +oo) by left-continuity. 

Equality (8.38) is proved as in Theorem 8.1, as well as the equality in (ev4). □ 

Remark 8.6. Using Remark 4.10 we deduce from the symmetry of that the solutions of 
the reduced e-regularized evolution problems with z^ly) > for every ye [— |, ^] satisfy 
the inequality zf{t,y) > for every t > and every y € [—5, 5]. Therefore any solution 
of the reduced problem for the approximable quasistatic evolution satisfies z^{t, y) > for 
every t > and every y G [—^,^]. In view of this property the hypotheses on Zq in 
Theorem 8.4 may seem artificial. This is not the case if, for instance, is replaced by the 
hexagon 

i^fce. :={(a,0GRxK:|a + CI<2, |a - C| < 2 , |5a-C|<7}, 
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with vertices (0,2), (-f,!), (0,-2), (|,-f). Indeed, let WT{t,y) 

2' 2J ' 



ipT{t)y for every t>0 and every y G [—5,5], where 



if < i < T , 
t iit>T. 



Let {ut, StiPti Zt) be an approximablc quasistatic evolution for the problem in dimension 
d — I corresponding to = ^ and K = Khex , with boundary datum wt and initial 
condition (uq, eo, po, zq) . We assume uq = 0, eo = 0, po = 0, and zo{y) = 60 for every 
y e [—5,5], where Co := — l^'(^o) € (~^)^)- By Remark 5.7 the functions exit), Prit): 
Zrit) do not depend on y . Wc can study the evolution of the corresponding fimctions crT{t) 
and Ct(0 Theorem 7.1 and plot them as paths of (arit), CT(t)) plane (cr, C) 

(see Figure 4). Since V'{6) = — 1 for 6 < —1, it is possible to prove that there exists T > 




Figure 4. A path (crT(i), CrW) from A := (0,Co) to P := (0,1) corre- 
sponding to an input (firit) , together with a path from A to Q := (0,-1) 
corresponding to a different input. 

such that {aT{T)XTiT)) = (-1,1)- It is then easy to see that (<TT(t), CtW) = (t-T-l,!) 
for T < t < T + 2 and (crT(t), CtC*)) = (1, 1) for t > T + 2 . Therefore any point at height 
C = 1 can be reached from any point (0, Co) with Co G through an appropriate 

loading- unloading path (see AP in Figure 4). Similarly, replacing (pT by —ifr, we can 
show that any point at height C = ~1 can be reached by another loading- unloading path 
(see AQ in Figure 4). 
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Note that, under the assumptions on zq considered in Theorem 8.4, the same oscihation 
phenomenon near y = occurs when K is replaced by Khex ■ 

Since Khex does not satisfy (6.5), the corresponding one-dimensional problem cannot be 
used to study evolutions of simple shears in elasto-plastic materials. However, they may 
be used in the study of uniaxial loading of cylindrical bodies, when the material exhibits a 
different behaviour in tension and compression (see Remark 6.2). 
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